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FOREWORD 


This la the final report; of a study made under Contract NAS 3-’2483 for 
NASA ~ Lewis Research Center, Cleveland, Ohio. Mr. Harold Schmidt at Lewis 
Research Center was technical monitor for the atudy. The report presents 
results of work performed during the 14 month period, August 1980 through 
September 1981. 

Ihe Lockheed -California Company was the prime contractor to NASA and the 
work was performed in the Commercial Advanced Design Division at Burbank, 
California. Important segments of the work which required special expertise 
were subcontracted to the Thermal Sciences Laboratory at Lockheed Missiles and 
Space Company, Inc., Palo Alto, California, and to the Visibility Research 
Center of the John Muir Institute, University of Nevada, Las Vegas, Nevada. 

The following individuals were principal contributors. 
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G. Daniel Brewer — Study Manager 
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John Muir Institute 
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ASSKSHMlvNT Of CRASH flRK 
HAZAUU Of Ul^ fUEI,El) AIHCRAfT 

0,1), Ri’ower, (1, WlttllU) li,f, Vcruaw 
K, Pfirmley*, K. Cir.ia*, f,(!. WoUho’t 


SUMMARY 

Liquid hydroKeii (LII2) la u in'lmfU’y ctmdJdntc fuel for future trtinsiuirt 
ulicraft. Studica performed for NASA liave bhowu Ita udvantapob in almoat 
every cate^'ory of comparison. AccardiuRly, when the international air trans- 
port industry is forced to turn to an alternative for conventional petroleum- 
base Jet A for reasons of either price or availability, or both, I.U2 is ex- 
pected to be a leading contender. 

In the meantime, there is much technology development which should pro- 
ceed in order that, when the time for a decision ie at hand, the decision 
regarding the choice of fuel for the Industry can be made on a credible basis, 
giving full consideration to all the factors which are important. One of the 
factors about which there is a serious lack of credible data is safety. How 
safe is liquid hydrogen compared to other fuels? What special precautions 
should be taken in design to accommodate it properly? In event of a crash 
in which passengers and crew survive the impact, would t*“*y be exposed to 
more or less hazard from fire or explosion if LH 2 is the fuel, rather than 
other fuels? - how about surrounding areas? What studies and experiments 
should be performed to resolve such issues? 

The subject study was performed to determine answers to these questions. 
Four scenarios describing situations involving fuel spills which can and do 
arise in commercial transport aircraft were studied to determine the hazard 
to airplane passengers and sur.oundings posed by LHg if it is employed as the 
aircraft fuel. The findings were then compared with results derived for liquid 
methane (LCH^), aviation grade kerosine (Jet A), and wide-cut gasoline (JP-4). 


* Lockheed Missile and Space Company, Inc. 
^John Muir Institute 
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In Humiiifiry I tm thp baHid oi tht uimlytilti ),l|^ vviii> hmiid to lu* 

tt HHfoi- lupl than any of tlip othtn* oand Jdatoo. Tlio paaaouKora on lu.ard. and 
l»t! 0 |>U« and pt’oporty In thu limiiodiate surronndlnHft, would both bo o/iitonod lo 
luHH hazard if u craabud niroraft; worn Inolod wUb 

Thin i'chuUh from tbrou fundamontal oonaldorat lonoi tin; 1,11^^ tanko nro 
leaa apt to auffor damnBo runultinR in npllU-d fiu-l j K npilUd, J,H^ vapor J, zoo. 
becomes buuyonit and dissipates into ulio aliiuniphort! so ivipldly U iiilnimlzoii 
the hazard; finally, 11 the vapor cloud Is Ifp'Ift-'d, tht- size and duration ol 
a fire from spilled are both siBulflcuntly less than with any otlier fuel. 

LH 2 tanks are considered leas apt to bu damaged because 

1) being mounted in the fuselage, LH„ tanks expose a far smaller 
dimension to frontal Impact, compared with conventional fuel which 
is stored in the wing and occupies nearly the entire span. 

2) they are designed to withstand a significantly lilgher pressure -tlian 
the rest of the fuselage and arc therefore leas apt to be the point 
of failure in a survivable crush. 

In the event a tank is ruptured and there is a large spill of fuel, 
does not spread as far as other fuels, it evaporates in much shorter time, 
becomes buoyant almost immediately, and finally, rises and dissipates into 
the atmosphere so rapidly that very little hazard is presented to surroundings 
outside the immediate area of the crash. For example, the following results 
were calculated for a survivable crash situation Involving spillage of nearly 

half the fuel load of aircraft designed to carry 400 passengers 10,190 km 
(5500 n.ml,) : 


Vapor Cloud 


Fuel 

Quantity 

Spilled 

Liquid 
Spread 
Rad lus 

Time 

to 

Vaporize 

Start 

to Rise 

Dissipated 

Time 

Distance 

T i me 

Distance 

llelj*hl 


(kg) 

(m) 

(a) 

(a) 

(m) 

(a) 

(m) 

(m) 

LII 2 

12,600 

35 

32 

12 

1.0 

146 

411 

575 


34,398 

61 

117 

1,560 

624 

1 ,624 

713 

72 

JP-4 

42,210 

143 

785 

Does 

not r:fs(’ 

11,736 

6,816 

0 

Jet A 

42,210 

331 

A 

! 



* 


*Thc vapor pressure of Jet A Is no 

u high 

enough at; 

standard atmospb 

it‘t'1 •' 

conditions to form a combustible 

mixture with the 

atmosphere. 
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I’liiany, in fivtint the spilled fuel la ignited, the duration of the F.ll^ 
fire would he so brief It would not Imat the fuselage to tlie point of eonapse, 
as would be the case with the otlier fuels, so the ]>nssengers would be safe If 
they stay In their seats until the fuel- fed fire Is burned out. I'or example, 

In the aircraft crash situation just cited, the resulting fire will burn Itself 
out In about 22 seconds In the case of MI2, vermis approximately A'l seconds 
for T.GII^, more than 10 minutes for , TP-4, and even longer for det A. 

The sub,1cct study represents a preliminary Investigation. The findings 
look very favorable for hydrogen. However, It must be recognized that there 
are many other cubuh and circumstances to be analyzed, the tools for analysis 
need to be improved, and experimental testing must be performed to establish 
a basis for validation of the computer model results before complete and final 
conclusions can be reached. 
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SYMHOI.a ANP UNITS 

» (louHCant 

2 

» uf iiul'f llriiugU wliicli Hlr in (sntvaliu'd (in ) 

■ uraa ol' liquid ixuil. (m^) 

■ ctmutunt 

« Hpeeiflc hcut 'iupacity ol; ulr at. conotaiU; pruHBuru (cul/gK) 

“ upuciric heat capacity of gaHuouu fuel ut conutant prcHHvive (cal/gK) 
■' specific heat capacity of puff: at constant pressuie (cal/gK) 

■ length parameter defining surface roughness (m) 

■ vapor pressure of water (mb) 

* saturation water vapor pressure In the puff (mb) 

« evaporation rate of liquid fuel (m/s) 

= Fanning friction factor for open channel flow 
■5 Froude number of the liquid flow 

■ acceleration of gravity (m/s ) 

“ enthalpy of air entrained during one time increment (cal) 

■ enthalpy of puff at start of a time increment (cal) 

■ latent heat of vaporization (kW sec /kg) 

■ length of puff in downwind direction (m) 

■ crosswind (cross runway) width of puff (m) 

■ latent heat of water vapor condensation (cal /kg) 

* latent heat of fusion of liquid water (cal/kg) 

■» maximum rate of vaporization from the liquid spread (kg/sec) 

- mass of air entrained during one time Inerement (g) 

® maximum imiss rate of vaporization per unit spill length (kg/scc*M) 
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nij, ■ nifiHH of fuel In puff (k) 

Wp " tiWHM of puff of HUut of time increment (jj) 

wmer entrained ilux'ing «)ne time Jnerumeut (ft) 

" iiioHH flowrute of uir (kK/fHoc) 
ni|P ■ mmo) flowrate of fuel va|>or (kg/Hoc) 

M " molecular weight ol fuel vapor or air 

" moleculur weight of air (g/mole) 

" moleculnr weight of fuel (g/niolu) 

■ moleculur weig.ht of puff mixture (g/mol.e) 

" moleculur weight of water (g/mole) 
n ■ Muimlugw roughiiunH coefficient 

f\ 

P ■ atmospheric pressure (dynes cm~^) 

|» “ exponent of height ratio in wind speed profile (-) 

3 

Qy ■■ total volume of the liquid spilled (m ) 

= volume rate of the liquid up ill (m /s) 

= volume of liquid spill per unit of spill centerline length (m'^/mi 

2 

q"^ “ convective heut flux from flames (kW/m ) 

q"p = incident radiation heat flux from flames (kW/m^) 

n 

q"j^ “ absorbed radiation heat flux from flames (kW/m‘‘) 

q" « heat rate per unit area to the solid surface (kW/m*") 

r »» puff radius (m) 

r^ B maximum radius of the spread from the center (m) 

R* B universal gas constant (erg/mole K) 

Re = Reynolds number of the flow 

t = time elapsed since start of vaporization from liquid (seconds) 
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tliiui to maximum axial aproad (a) 
time Incremant (a) 
compernture of air <K) 
flame temperature (K) 

puff temperature at a tart of time Increment (K) 

Intermediate puff temperature after balanclnR aenalble heat 
tranafer (K) 

Intermediate puff temperature after condenaing water vapor (K) 

intermediate puff temperature after freezing water dropleta (K) 

Intermediate puff temperature after any melting of Ice cryatala 
and re~evaporatlon of water dropleta (K) 

wind apeed (m/s) 

horizontal downwind apeed of puff (m/s) 

entrainment speed (m/s) 

first term of entrainment speed (m/s) 

second term of entrainment speed (m/s) 

aircraft velocity along spill line (m/s) 

volumetric flowrate of air at atmospheric temperature and 
pressure (m^/s) 

mean horizontal velocity of the liquid flow (m/s) 
volume of puff (m ) 

volume of puff at end of time increment (m^) 

volumetric flowrate of fuel at atmospheric temperature and 
pressure (m^/s) 

Initial puff volume (m^) 

volume of air entrained during one time interval (m^) 
rate of fall of still liquid surface (m/s) 
vertical puff speed (ra/s) 



• distance from the center of the spill line (m) 

■ maximum lialf axial spread • from spill centerline (m) 

■ downwind distance of puff from origin- (m) 

" change In downwind displacement (m) 

- minimum combustible mole fraction of fuel in the uir 

• height of puff above ground (m) 

- reference height for wind speed profile (m) 

- height of puff above ground (m) 

■ puff thickness In vertical direction (m) 

■ lapse rate of temperature in the atmosphere (K/m) 

• liquid flow depth (m) 

■ effective flame emissivity 

- density of air at height of puff (g/m^) 

■ density of gaseous fuel at 1 atmosphere pressure and boiling 
point (g/mJ) 

■ density of liquid fuel (g/m^) 

■ density of puff (g/m^) 

■ wave velocity (ra/s) 

■ liquid viscosity (kg/s m) 


SECTION 1 


INTRODUCTION 


This la the final report of a study performed for NASA-Lewls Research 
Center to assess the hazard posed by liquid hydrogen when used as fuel in 
conanercial transport aircraft. In addition to liquid hydrogen (LH 2 ) , the 
hazards presented by three other fuels were also calculated for comparison. 

The other fuels were liquid methane (LCH^) , Jet A (aviation grade kerosine 
currently used in commercial transport aircraft), and JP-4 (the widecut 
gasoline used as jet fuel by the Air Force). 

Four specific scenarios or descriptions of incidents leading to spilled 
fuel were specified by NASA to serve as a basis for the safety analysis. The 
incidents ranged from a very small leak such as might result from a crack in 
a faulty weld on a fuel line or tank, to a catastrophic crash which would 
cause simultaneous rupture of all the fuel tanks on board the aircraft. 

The study was initiated because is a candidate fuel to be used when 
Jet A, produced from crude oil, is no longer available to the world transport 
fleets on an economically acceptable basis. Previous studies conducted for 
NASA have shown numerous advantages which can be gained by using 
These include 

• Global availability . It can be made from a renewable resource, water, 
using any source of energy which might be locally available. Fossil 
resources are not required. 

• Pollution. Its major product of combustion is water vapor. It emits 
none of the noxious effluents of carbon-containing fuels save a mini- 
mum amount of NO , 

X 

• Noi^. Because they are lighter and therefore can use smaller engines, 
I.H 2 "fueled aircraft would be about half as objectionable to airport 

ne ighbors. 


8 


• «ii tixnmplo of an energy resource which can be used 
to produce any of the alternate fuels . Fewer tons of coal would be 
required to produce the fuel used by aircraft to fly their design 
miss ions if MI2 is the fuel. 

• Direct Operating Coat. hH2''fueled aircraft can bo conipetitlve in DOC 
using current fuel production technology^ ./ 1 th advanced processes • 
I.H2 can be produced from coal and water at costa which offer signifi- 
cant advantage. 

The relative safety of LH2 has been subject to question. The present 
study Is an important step toward establishing a sound basis for judgement 
on this matter. In addition to providing a preliminary answer to the question 
of how LH2 compares with the other fuels in terms of crash safetyi it also 
provides suggestions for additional analyses and experiments which are con- 
sidered necessary for final resolution of the problem. 
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SUCTION 2 


TECHNICAL Ain'ROAClI 


2.1 Study Team Organization and Work Plan 

The study team consisted of representatives from three companies selected 
to provide expertise in the diverse technology .ireas needed to perform this 
study effectively and efficiently. 

The Lockheed -California Company was the prime contractor to NASA and 
managed the study. LCC also reviewed transport aircraft accident data from 
various sources to derive a basis for establishing what would constitute 
credible damage associated with each of the crash scenarios specified by NASA. 
These damage assessments were then considered in conjunction with aircraft 
designed for each of the subject fuels to specify typical fuel spill rates 
and quantities for each scenario. 

Under subcontract to LCC the Thermal Sciences Laboratory of Lockheed 
Missiles and Space Company, Inc. (LMSC) was responsible for analyzing the 
hazard resulting from spill of the liquid fv:;.LS. A computer model was created 
by Mr. Richard Clma to characterize the spreading and vaporization of volatile 
liquid spilled in both a radial and an axial pattern. The latter represented 
the spill situation where, for example, an aircraft strikes an obstacle during 
a landing shortly after touchdown and spills fuel during deceleration to a 
stop, A radial spill pattern would be produced when impact with an obstacle 
stops the aircraft Immediately. 

The geometry and mixing characteristics of the cloud of vaporized fuel 
as it disperses into the atmosphere were modeled l)y Dr. lirlc Walther and co- 
workers at John Muir Institute. The geometry of the cloud, the profile of 
the flammable mixture within it, and the path followed by the dispersing cloud 
under various atmospheric conditions were needed to provide information about 
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the extent and duration of the expoauro to hazard of the aircraft pnaHengerB 
and the area aur rounding the craah acene, for example, the tlme^hl story of 
the geometry and location of the flammable mixture within the dispersing cloud 
provided Information about the potential size of fire which might result If 
the cloud were Ignited at various clmea. 

Another type of hazard which was explored by LMSC related to the problem 
of how to handle a small leak of fuel into an internal cavity in the airplane • 
There Is the possibility that if such a leak occurs, an explosive mixture 
could develop and a detonation result. The problem was resolved by determin- 
ing, for various size leaks, the quantity of air flow required to assure that 
a flammable mixture could not develop in the compartment as a whole. A 
flammable mixture will always exist locally at the leak source of course. 

Finally, hazard to passengers inside the aircraft cabin was evaluated for 
the case where the outside of the fuselage is exposed to flames from burning 
fuel spilled as a result of one of the survlvable crash scenarios. 


2.2 Study Guidelines 
The objectives of the study were: 

• To assess the hazards associated with non-normal and crash landings, 
failed takeoffs, and ground accidents ln'''olving turbine-powered air- 
craft fueled with liquid hydrogen (LH 2 )> 

• To compare LH 2 -fueled aircraft with similar designs fueled with liquid 
methane, Jet A, and JP-4, on the basis of: 

- Post-crash fuel system damage 

- Potential deflagration and detonation hazard 

- Survivability potential for persons In and adjacent to the damaged 
aircraft 

• To formulate a program of research and testing found to be needed to 
fill technological gaps. 

The hazards sssoclated with LH 2 were to be evaluated end compared with 
those from liquid methane (LCH^), Jet A, and JP-4 in connection with the 
following accident scenarios: 
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1. A non-normnl landlnR or Rround necldent wliioh results in fuel system 
insulation danuiRO and/or fuel system darnsRo permlttinR liquid hydro- 
Ren to vent, eseape, leak or run out of a punctured tank or broken 
line. 

2. A survivablc crash landlUR or failed takeoff where damoRe to fuel. 
tmikuRe or lines results in massive release of liquid hydrogen after 
the aircraft has come to rest. 

3. A survivablc crash landing or failed takeoff where damage to fuel 
tankage or lines results In massive release of liquid hydrogen upon 
Impact and during aircraft deceleration. 

4. A nonsurvlvablc or catastrophic crash resulting In the maximum rate 
of energy release in the form of a conflagration and/or explosion. 
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SliCTlON 3 


INPUT DATA 


In thia acction, data rcquirud ua input to perform the crouh haisard 
aimlyaea of eadi of the four fuela in the altuationa specified by NASA are 
presetiLed. In turn, aubaections present data on pertinent properties of each 
of the fuels; results are described of previous experiments with LH^ and LCU^ 
which were performed to assess some of their safety-related characteristics; 
typical designs of transport aircraft for each of the fuels are described 
which were used as a basis for establishing damage likely to result from the 
specified crashes; results of study of aircraft accident records are reported; 
and finally, the damage selected as representative of that likely to occur to 
the specified aircraft in the subject crash scenarios is described, along 
with resulting fuel spill rates and quantities. These fuel spill data are 
used in work described in subsequent sections in evaluation of the consequences 
of the spills. 


3.1 Fuel Characteristics 

A number of properties of the four fuels considered during this study 
are presented in table 1. The properties selected for presentation include 
some which are of general interest in comparing fuels for the aircraft appli- 
cation; most are properties which were involved in the analysis of hazards 
associated with spills resulting from the kind of aircraft crashes postulated 
for this study. 

As noted on the table, the properties of the cryogenic fuels, and 
LCH^, wen: t; ken from reference 5, a recent compilation by the Center for 
Chemical Engineering at the National Bureau of Standards in Boulder, Colorado. 
The data for Jet A and JP-4 were primarily derived from the respective speci- 
fications for the two fuels, references 6 and 7, except as noted in the table. 
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TAULK 1 


siiU-;cTKi) pRoj'jiiiTiiis or hikls 



Hydrogen 

Metliuna 

Jet A 

JP'4 *^* 

Nominal Cumimtitlon 

•‘2 

CII4 

‘^'*1.93 

0'*1.93 

Molecular Weight 

2.016 

10.04 

- 108<®> 

-132*0* 

Heat of Combuttion (low), h J/g 

120 

50.0 

42.8 

42.8 

Liquid Ueniity, g/cm^ at 283 K 

0.071* 

0.423* 

~ 0.811 

-0.774 

Bolling Point, K (at 1 Atmoiphoro) 

20.27 

112 

440 to 539 

333 to 519 

Freezing Point, K 

14.4 

91 

233 

215 

Specific Heat*, J/g K 

9.69 

3.50 

1.98*®^ 

2,04*0* 

Heat of Vaporization, J/g 
(at 1 Atmosphere) 

446 

510 

3G0*®* 

344*0* 

Vapor Pressure, kPa 

- 

- 

- 

18 

Diffusion Val. in NTP Air, cm/s 

'::2.00 

<0.51 

<0.17 

<0.17 

Buoyant Vel. in NTP Air, m/s 

1.2 to 9 

0.8 to 6 

Non Buoyant 

Non Buoyant 

Flaminabllity Limits in Air, Vol. % 

4.0 to 75.0 

5.3 to 15.0 

0.6 to 4.7*®* 

0.8 to 5.8*0* 

Vaporization Rate W/0 Burning, 
cm/min (From Liquid Pool) 

2.5 to 5.0 

0.05 to 0.5 

N/A 

0.005 to 0.02 

Min. Ignition Energy in Air, mJ 

0.02 

0.29 

0.25<*“> 

0.25***** 

Autoignition Temp., K 

858 

813 

>500 

>500 

Rate of Liquid Lowering of 
Burning Pool, cm/min 

3 to 6.6 

0.3 to 1.2 

0.17*®* 

0,20*0* 

Burning Vel. in NTP Air, cm/s 

265 to 325 

37 to 45 

18**0* 

381**01 

Flame Temp in Air (Stoichiometric), K 

2318 

2148 

2200**0* 

2200**0* 

Thermal Energy Radiated to 
Surroundings % 

17 to 25 

23 to 33 

30 to 42 

30 to 42 

Detonability Liniits in Air, Vol. % 

18.3 to 59.0 

6.3 to 13.5 

N/A 

1.1 to 3.3 


"At Normal Boiling Point 
NTP - Normal Tern pirature and Prestura 


3.2 PrevlquH Satety KxperliwuntH 

Thqrq ia a flurin’isingly inq«B«r lilst^ *y of wuIl^-oigmiiKod, sclootlf ioully 
uxucutod safoty oxiierfiiioiiCfl conduetod with elthor Ul^ or to ootablluh 
the UHture uud extent of Imisnrdo nnHoclnted with large iiintantmioouo HpiJlH, 
or long dnrntion spiUh of nlxablo quautltien. 

A aeareh waa mode tor experimuntul cvtdenee to be used for validation of 
the unulytleaX modela eruatod for uae in the aubjeet atudy. in the limited 
time available, the following four aourcea were found to provide the moat 
uaeful information: 

• Teatlng performed by the Advanced Development I'rojecta Organization 
at LocWieed -California Company in 1956-57(11) . 

This very early work at Lockheed's "Skunk Works" included a series of 
teats to determine the detonability of LII 2 . Although the tests involved only 
small quantities of LH 2 . they were informative in establishing that LH, would 
not detonate, even when solid oxygen was Immersed in the liquid hydrogen, 
unless a strong explosive charge was used as an initiator. 

In the series of experiments Lil 2 contained in a thermos bottle was sub- 
jected to heavy impact 61 times. Ignition never occurred as a result of the 
Impact alone. Wlien a squib or a hot wire was employed to ignite the hydrogen 
vapor following smashing of the thermos bottle by a heavy weight, very little 
overpressure was measured, indicating deflagration without detonation. De- 
tonation occurred in two tests, both in cases whero oxygen was deliberately 
mixed with the LH 2 . 

• Spill tests performed by Arthur D. Little, Inc. for the Air Force 
in 1958(12). 

This teat program was designed to simulate a set of conditions whlcl: might be 
encountered In storage or transport of large quantities of LH 2 , with the ob- 
jective of establishing a basis for realistic quantity-distance relationships 
and safe-handling procedures. The tests included spills of several quantities 
of LII 2 ranging in size from 5 liters to 5000 gallons to study the magnitude 
of the hazards involved. Spill tests of some hydrocarbon fuels were also 
made for comparison purposes. Once again it was demonstrated that hydrogen 
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will uol; Ignite from Impaet, either from bulletH being fired through the 
contolner or from the container being Hinaabed, and thut. uiiIohn confined on 
at leaut three nldea, it will not detonate when Ignited, 


The comparlHon tentn with hydrocarbon fnelii demonatrated the difference 
in burning chanicterlatlcH, Spills of 32 gallons of each of the fuels pro- 
duced the following resultH: 


Fuel 


Duration of Fire 


LH2 

Propane 

Gasoline 

JP~4 


27 sucotids 
4 minutes 
^ 5 minutes 
minutes 


Radiation measurements of the hydrogen flame indicated a maximum emisslvlty 
of 0.083| compared with a value of about 1.0 for the hydrocarbon fuels. 

Following are conclusions taken from reference 12 which are considered 
pertinent to the subject study; 

1, Detonation effects of a large-scale spill - There was no evidence 
to substantiate the premise that the mixture of hydrogen and air 
formed by the sudden release of liquid hydrogen will detonate when 
ignited by a spark. 

a. In no spill did the resulting vapor cloud detonate. 

b. In all tests (In which an ignition source was provided), ttie 
pressure effects of the resulting fireball were negligible 
except where there was Beml-confinement of the gases. 

c. Transition from deflagration to detonation requires a significant 
pressure buildup. 

d. Detonation of a perfect stoichiometric mixture of hydrogen and 
air took place only when Initiated by a strong detonator. 

e. It Is probable that a perfect lulxture of hydrogen and air will 
not occur. 
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2. ))t> n /»gra|. loii^ t' f fc<t.i K of I argo-MPolo 1 J<iu .ld"hydr oRon rp.UoiHP - Xlip 

Fas'.Hrda atiHopimoti with t he oj’Pun Jng aTlur tlu* j.ol oatoJ iiY 1 lij» J d 
hydrogou voiy /u.pordlug ti* tho oiul dttroiJoii of tho (lainc-n. 

Oi lUd J agral ion oiit-'otn oro oonHidoralily loan li Ignition oppuj'm 

liimiodintoly niion l linnUI riflooHOf initial l iroboll niiu' Im'vo/mtsi 
grcM' ly with doloyiul Ignit ion. 

h. in (ill invgi' npitih o) hydt’ogv'n, tlio din’otiini oi tlio l irohaJi 
wuuid ho no hrUd that; inotal atria'lnroH would mil l oi; no nlgnili- 
I'anl datiiugo ; howovoi'i coiiihniit Ihlo itiiitorlaJn would probably bo 
Ignittid and pomonnol within tlio liroball would bo n.^oroly 
buraod. 

c, Al ter diDsipat ioii oi the Initial fltoball, I lamon would bo coil'’ 
iiiiod to a region directly over the diked area, but could extend 
150 ft into the air. 

Hydrogen haaar dH an com|aired to thoae of liydrocurbonH " The hazard a 
anboclated with the handling and atorage of liquid hydrogen are leaH 
than thoue encountered with petrolenri fuela, 

a. Liquid hydrogen evaporateu much mure rapidly. (When a flame waa 
present, a pool of iiydrogeu vapor ifted in 1/20 to 1/50 the time 
required to evaporate a gaHollnu or propane pool of equal depth.) 

b. The flame from combustion of liquid hydrogen radiates less than 
1/10 the energy per unit flame urea emitted by the hydrocarbon 
fuels (cmisslvity of the hydrogen <0.1; emisslvity of a hydro- 
carbon vl) , 

c. The high evaporation rate and low heat transfer from flame to 
liquid pool make che probability of boil-over (i.e., the super- 
heating of a large mass of liquid until a sudden flashing to 
vapor by a sizable portion of the mass occurs) Imi^robabLc with 
liquid hydrogen. 

4. Effects of storage-tank rupture - In cases where rupture of u storage 
tank Is caused by a pressure buildup in excess of the capacity of 
the pressure-relieving device, a fire may be caused by the tearing 
of the metal at the time of rupture. 

• Spill tests of LU., performed by NASA at White Sands Test Facility, 

New Mexico, during 1S'80(^^2, 

3 

These consisted of seven spills, six approximately 5.7 m (1500 gallons) of 

3 

LH 2 , the seventh was about 2.7 m . No ignition was planned and none occurred. 

The objectives were to determine the time required for evaporation of the LM^ 

and to atudy the path of the gaseous hydrogen as it dissipated In the 
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auiioHphpru. l»j pi linlnftry duni, hup|ilftmontpd Uy moviu atul will I p'nni up,ra|>hfi 
hlutwliif.; iliP movomPiu ol t tiP clmul oi watnj; vii|iur Inniifi! . 11 ; Ou' tvu.i'tjui. Iivdrn^'pu 
dittaipjiiud, iirovldpd n banif* fpr v^Uddi luK iIip mailipmat U al ihoiIpI uf auch an 
aviait di>vplo|u'd J <m’ t lie aulijiTi atudy. 

• il!)Jll o( ll(|uJd nalnral «'oiulmapd by l.awrciai' I, I viTtiioH' 

National Laboratory l or Hopartiiiont ol Km’rj’.y al (IbJita t.abo . 

Luwrmico Livoniioro National L/iboratitry (I.LNL) la fombn t lay, an oxloiialvo 
Huritsa ol.' Hplll toHl.a with l.N(; uiidor oontraot to tho Uopart men; ol r.ooi yy 
(UOU), Thu piupoau of thu work la to dt!Lc:rmlm! oKpur Iniontal 1 y bow LN(J dlaal- 
patoM in the atmobplujrei und tbui to rupruHont thu ocuurroiK'.o with a niatlu" 
tnutlcal model. The objoctivu iw to provide DOli with tho oapablllty ol prudio.t" 
ing the huzurd to burrouiidiiig arouy in the event large quuntltlea ol l.Nti are 
apiU.ud. Since natural, guu la t-jenerally at leant 85 pereent methane the re- 
aultb of LLNL’s teat program were of inturetit in determining the validity of 
the atmospheric dispersion computer model created by .lohn Muir Inistitute lor 
the subject program. 

Two series of LNG spills hove been conducted by LLNL at China Lake, 

3 

California; an initial series used 5 m in each spill, a later series in- 

3 

volved nine spills of approximately 28 m (rouglily 7400 gal.) each, The 
tendency of the visible cloud associated with vaporized LNC to remain non- 
buoyant for extended periods compared to the relatively rapid rise oi the 
cloud resulting from a LH 2 spill as observed in the NASA experiments provided 
confirmation of data calculated by thu subject atmospheric dispersion model. 

These four experimental programs were helpful In providing information 
which was of use in carrying out the evaluations in this study. However, it 
is obvious much more experimental testing should be done to provide authentleu" 
tlon of thu analytical conclusions presented herein. Section 7 of this report 
suggests o number of examples. 

3.3 Aircraft Design Characteristics 

During previous studies for NASA(^), I-ockheed established preferred de- 
signs of aircraft fueled with LII 2 , LCH^, and .let A, respectively, for a series 
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gf Alrt’-raJU dtuilf?mul Cor oiig gt tliiiHc* iiiljui loiiw wen' (ipU'ptgd to 

Hgrvg «« <i buHig lor Mu? cv/ttiU fin* liaj^ard cniniair I huh oI tlig [ivtaipail nltuly. 

Tlip iniHMluu Hulypt-gil wan tu parry /i00 |tan(.(>nH«'rh i plan ap|<ru|U lau* carKa, a 
UualHu rnup,g o< 10,190 km (-WOO ai a nominal criUm* Hiaaul iil M/u'li O.H't, 

tinmu t'liar/icter liil Inn ol tla^ aJrnran dan ly, nod Idr i'aoh n1 llip npnc IJ li d 

iup>,la to acu'.om})!l .i.idi t.lnia mlaaJon arn ahown In tab In II. lUu.auiio ni ilin a.lm.1-" 

l.arlty of the huatinp, valuiu) and pliyoilnal propnrUoa al Jot A and JJ'' A Jt wan 
asHumud, for puritoacn of tbla ntudy, that Llin diaraptorliitlpn lintod Jar tin* 

Jut A-fuulud aircraft al.nu rupram-'iit tbu Jl’-Jl du)i,ly,ti. 

TAULD 2. ” CIlAUACTblUaTJCS OF Ul„, I.Clt., JKT A, AND JP-4 FUKl-FD 
SUUSONIC TUANSl'OimJ ^ 

(400 Fanaununr, 10,190 km (5500 n.tni.), Mach 0.H5) 





Liquid Hytlrugun 

Liquid Muthuiio 

Jut A und JP 4 

Grou Wsight 

kg 

(lb) 

1U8,740 

(372,000) 

22b,680 

(497,300) 

232,060 

(611,0110) 

Tutal f uul Walglit 

ko 

(lb) 

2li,6U0 

(66,460) 

69,040 

(162,200) 

84,700 

(166,900) 

Block Fuel Weight 

kg 

(lb) 

21.G20 

(47,670) 

58,960 

(130,030) 

72,360 

(169,500) 

Oiieratlng Em|)ty Weight 

kg 

(lb) 

103,300 

(227,760) 

116,170 

(266,120) 

107,360 

(230,700) 

Wing Area 

2 

m 

(ft*) 

296.7 

(3,196) 

386.0 

(4,144) 

3B0.2 

(4,093) 

Span 

in 

(ft) 

61.8 

(170.0) 

6G.9 

(103.1) 

60.6 

(182.0) 

Futelege Length 

in 

(ft) 

66.7 

(216.6) 

61.4 

(201.3) 

60.0 

(197.0) 

Fuielage Diameter 

m 

(ft) 

6.63 

(21.76) 

6.10 

(20.0) 

6.84 

(19.17) 

Thrust Per Engine 

N 

(lb) 

136,000 

(30,350) 

177,030 

(39.800) 

106,030 

(41,600) 

Specific Fuel 

kg/hr 

“IT 

Ib/hr 






1 

Cuntumption (Cruise) 

lb 

0.0206 

0.202 

0.0604 

0.494 

0.0616 

0.003 1 

L/D (Cruise) 


.... 


17.4 

19.21 

19.13 

Takeoff Field Length 

in 

ft 

2440 

8000 

2430 

7973 

2431 

7070 

Approach Speed 

m/sec 

kts 

71.0 

138 

66.6 

129 

66.3 

127 

Deta from Reference 4 






Tlio p.ooimary mid ponural. arranpi-Mnmit nf the throe atrcraCt dealgne arc 
lUuHtratod in tlputoH 1, 2 , mid 3. It will he noted that the tanka for the 
eryopevilually'^’fueled aircraft: are Jocated wltliln the fuaelape, one junt aft 
of the f I Ipht Htation and oiv. Juat forward of the empennaEO. An extenalve 
detilpn Htudy was made for both the Ul2 and the hCH^-fuelcd aircraft to deter- 
mine If other tankage arrangomontu might offer advantage. It was found that 
I'or both fuels the fore and aft internal fuaelage arrangement provided the 
beat arrangement; i.e., lowest direct operating coat for a practical design. 
Accordingly, the tankage arrungementB shown in figures 1 and 2 were adopted 
for the LII^ and LCH^ aircraft designs, respectively. 

The noncryogenic fuels, viz.. Jet A and JP~4, are conventionally carried 
as snown in figure 3; i.e., with the fuel contained within the wing structural 
box. 


3.3.1 Fuel system arrangements 

3. 3. 1.1 Conventional fuel system. - Figure 4 is an illustration of a fuel 

system arrangement for a conventionally-fueled modern transport aircraft. 
Insofar as possible all fuel system plumbing Is carried within the wing struc- 
tural box in the same volume in which the fuel is contained. Thus, in event 

of a leak in the plund)lug no fuel Is lost. Also, burying the plumbing within 

the tank eliminates concern about damage to that part of the fuel plumbing 
system in event of a crash. The portions of the plumbing system Inside the 
wing box can be damaged only if the tank itself is penetrated, an event of 
much more serious consequence. The only portions of the fuel plumbing system 

which are vulnerable to primary damage leading to large fuel spills In event 

of a survlvable crash are those lines and components external to the tanks. 

For example, the line leading from the tank to an engine will be sheared off 
if the engine pod is torn from the wing. 

In a four-engine aircraft the auxiliary power unit (APU) is usually 
mounted in the tail section. Accordingly, a fuel line must run from the wing 
tank aft through the fuselage, back to where the APU Is located. Under normal 
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fueled traimpurt 







Figure 4. - Fuel system arrangement for conventionally fueled transport aircraft. 




circumatancc's this fuel line cun be located well within the fuselage structure 
ami thus bo virtually invulnerable to serioua damage, except in cases where 
the fuselage is physically torn apart. For example, in the three-engine 
l.ockhued L-1011 aircraft, the fuel line which runs aft to the No. 2 engine 
and to the APU, is located Just below the cabin floor and very close to the 
vertical centerline. There is abwmt eight feet of structure below Che fuel 
lino to the aircraft skin which effectively protects the fuel plumbing In the 
event the airplane makes a wheels-up landing. 

In addition to being located in the fuselage where the line is physically 
protected, since 1967 fuel lines in commercial transport aircraft have been 
required by Federal Aircraft Regulation to provide "a reasonable degree of 
deformation and stretching without leakage" to accommodate limited amounts of 
fuselage separation. Again as an example, Lockheed policy is to make the 
fuselage-mounted fuel lines capable of stretching to 50 percent of their 
original length. Thus, it would require a crash severe enough to separate 
sections of the fuselage by several feet Jn order for the No. 2 engine/APU 
fuel feed line in the L-1011 to pull apart and spill fuel. The accident data 
analyzed and reported in section 3.4 shows the probability of this happening 
in a survivable crash is fairly small. 

3. 3.1. 2 Cryogenic fuel system. - The general arrangement of a fuel system for 
a cryogenlcally-fueled altplane Is depicted in figure 1. With fuel tanks 
located fore and aft in the fuselage it is apparent fuel lines must run be- 
tween the tanks and also out to the engines mounted in the wlngn. In this 
case, as contrasted with the conventionally-fueled aircraft, even though the 
lines running to the engines are located In the wing structural box for maxi- 
mum protection against damage, they will be the source of fuel leakage result- 
ing from any Impact that severs or massively punctures the wing. 

The comparison of fuel spillage which might result from wing damage In 
a cryogenlcally fueled airplane versus one fueled conventionally is a matter 
of quantity and fuel spill rate. With Jet A or JP-4 fuel, both of which are 
stored in the wing, a puncture of the wing box will result In spill of all 
the fuel In that tank at a rate which Is a function of the size of the hole. 
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With ov LCU^, thu worat that will liupinm In uvmit; of wIuk tliima};i.' in tlmt. 
un onginu food lino will bo aovorod. Tlio aimmiu, of fuo.l apl,l ,U>d an a iroriul t 
will be contingent on whotber tlui booat pumps and t.lin fuel ahut’off valve are 
tunied off ov not. The rate of spillage will almoSL certainly bo. far loss 
than in the case of the conventional fuels. 

fuel lines for cryogenic fuels will be insulated to niinimliie heat leak. 
Without adequate insulation, the fuel would vuporiiio on its way to the engine 
fuel pump, ultimately causing the pump to cavitate, thereby starving the 
engine. Accordingly, the lines will typically be designed as shown in fig- 
ure 5. The line carrying fuel to each engine will be approximately 2.54 cm 
(1.0 inch) in diameter for the size aircraft adopted as a basis for this 
study. It will be encased in a closed-cell, rigid plastic foam about 3.81 cm 
(1.5 inch) thick. In turn, the foam insulation is contained within a thin 
aluminum tube 10.16 cm (4 inch) in diameter which serves the dual function of 
mechanically protecting the foam and preventing air from cryopumping In event 
of local failure of the foam. Bellows will be provided at intervals to ac- 
commodate thermal expansion and contraction, plus mechanical vibration and 
flexing. Wherever these fuel lines pass through pressurized areas they are 
contained in a shroud which is provided with positive air flow. Hydrogen 
(or methane) detectors will be located at the discharge to sense leaks. 

For a more complete description of a representative cryogenic fuel sys- 
tem see reference 3. There are many more elements to a complete system which 
are required to perform the various functions needed on board an aircraft. 
These include pressurization, venting, fueling, defueling, crossfeed, etc. 
However, for purposes of assessing the crash fire hazard, the engine fuel 
supply system just described, represents the principal challenge. Crash 
damage to the other systems, with the exception of crossfeed, will not result 
in significant fuel spillage. 

The crossfeed system involves the fuel lines which traverse the fuselage 
connecting the tanks so that there are multiple routes by which any engine 
can be fed fuel from any tank. They are located in the fuselage Just below 
the cabin lower floor as shown in the right side of section AA in figure 1. 
These lines arc subject to damage in event of a crash in which the fuselage 
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1 



Figure 5i - LH2 fuel line cross section. 

is broken. They are the same size as the engine fuel feed lines in the wings 
and their severance will result in the same fuel spillage potential. 

Cryogenic -tankage. - Tanks for cryogenic fuel are designed to entirely 
different specifications than their counterparts for conventional fuels. Both 
LH^ and LCH^ are tanked as saturated liquids at a nominal pressure of 145 kPa 
(21 psi). The tanks are pressurized with their own vapor product, and so they 
contain no air or oxygen. Consequently, there is no fire or explosion hazard 
v;lthin the tanks. Conventional fuel tanks, on the other hand, are vented to 
the atmosphere to minimize differential pressure on the structure, so there 
is continually a mix of fuel vapor and air above the surface of the liquid 
fuel which may be susceptible to ignition. The cryogenic fuels therefore are 
inherently safer in this regard in event of a crash. 

Tank pressure in the case of the cryogenic fuels Is controlled by venting 
to relieve excess pressure when heat buildup Is too great, and by intentionally 
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vapoviulu^ lliiulil liu’] wht'.n tank priuiHurts btiluw a appcified limit. TIkj 

cryof-itmic tankH arc duHi^iu’d to mi ultinmto iu’ohhui’o of 208 kl’u (80.2 i>h1) . 
lor comparison, rlu> pasHciifior cabin Is deu4'nod to an ultiniato proHsuro of 
only 12/* kPa (18 psi). n. ,i.a much more Ukoly, thoroforo, that the fuel tanka 
would Hurvlvo intact in ovont of a crash, rather than the rust of the fusolago. 
The fusclaBo would tend to break In its weakont parts. 


Another aspect of the tankage of the cryogenic fuels which results in 
tlieir being less likely to auffer damage in event of a survlvable crash, com- 
pared to conventional fuels, is tlieir location in the fuselage. The tank di- 
mension (its width) which is vulnerable in a collision with an obstacle on the 
ground is significantly less for the fuselage-mounted tanks than for wing tanks. 
As an example, in the LH2“fuelcd aircraft shown in figure 1, the fuselage tank 
has a width of 6.63 m (21.75 ft). The Jet A-fueled aircraft designed for the 
same mission, illustrated in figure 3, has a wing span of 58.5 m (192 ft), 
most of which contains fuel and is vulnerable to damage from impact which 
could cause fuel spillage. In the LH^ aircraft, of course, the wings are 
devoid of fuel except for the engine supply lines which, as explained in 
section 3. 3. 1.2, are contained in the structural box between the front and 
rear spars and are downstream of shut-off valves controlled from the flight 
station. 


An additional factor in this regard is the consideration that the fuselage- 
mounted tanks of the cryogenic fuels have a significant amount of structure 
both ahead of and beneath them to absorb impact loads. Examination of fig- 
ures 1 and 2 shows that the nose structure of both of these aircraft would 
have to be crushed before a frontal collision would impact the forward tank. 
Also, the underside of the tanks in both aircraft is protected, by a minimum 
of 45.7 cm (18 in.) for the forward tank and 35.5 cm (14 in.) for the rear 
tank, with specially designed structure to preserve their integrity in event 
of nose gear collapse, a wheel s~up landing, or a tall scrape. This protection 
should assure a better record of crash survival for cryogenic fuel tanks than 
the conventional fuel wing tanks. 

The tanks for both LIl^ and LCH^ are insulated to minimize heat leak 
which would cause rapid boiloff of the fuel. Two types of insulation which 
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Imvu Utiui louiid to l»o moot promlolaB h» « jccault of dutuilod 
utti illuwtvatud in fiBovvi 6 oud 1 oh doHlyuod fot tliu Lllj^-luolod olircraft. 
ThoHO luHul atiou uyattiiim art*, rnapccr Lvo ly , rigid cloHod-iinl I J cmm, and 
mlci'otipluo’nn. 

Both doHlgm! wort* nnUibl. lalit'd on tho banlw of wlii:hiil?.lnB dlrt'ot oj)i*ra( lug 
t'.OHt for tlio tili'oraft by holcctlng an t)jjtlinum comproml ho botwoou iuHulatlmi 
KyHUeiii wulRlU, and coHt of fuol. nllowtid to boll, oft. Tlio rigid cloood-oell 
foam dooigu muot bo tliJckor Limn the microupboro HyKtom botmuao. the ovuouuLod 
inlcrouplu'ro aunuluo in a more cfflolont Inoulutor. Uowovur, tlio. mlcruupluiro 
sytitem weigha a little more and ruquirua an active pumping ayatem, an extra 
complication not rcquiretl by the paaaive foam aystem. 

NASA haa carried out extensive preliminary tests of many candidate foam 
InHuIauiona''^^. Some, such as Stepan f’oam and a General Electric polyurethane 
composition, have been shown to bo very efficient for the cryogenic applica- 
tion. The potential life of these materials has yet to be explored, however. 

As shown in figure 6 the rigid closed-cell foam is applied directly to 
the outer surface of the tank. Over that is laid a vapor barrier such as 
a 0.127 mm (0.005 in.) thick composite sandwich of mylar-alumlnum-alumlnum- 
mylar-fabric mesh (MAAMF) , the purpose of which is to prevent cryopumping of 
air into the foam during its operational life, and also to prevent the same 
occurrence in the event the fodm develops a crack or other defect. 

A flexible open-cell foam layer is applied over the MAAMF to provide a 
spring-like cushion between the rigid foam and the Kevlar-syntactic foam panel 
which serves as the outer surface of the aircraft In the area of the cryogenic 
fuel tanks. A second MAAMF layer le applied between the flexible foam cushion 
end the Kevlar panel to provide a sealed annulus occupied by Che cushion which 
is purged with nitrogen. 

The microspheru insulation system is illustrated in figure 7. Micro- 
spheres are tiny hollow spheres of borosilicatc glass. They are contained in 
an evacuated annulus surrounding the tank which is formed by a 0,127 mm 
(0.005 in.) thick stainless steel (ORES) outer shell. The annulus is 
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Figui'c 6. ~ Rigid, cloeed-ctill plastic foam tank insulation 



Figure 7. - Evacuated microsphere tank insulation system. 



nmlntalnod at n aqft vacuum of approx I matoly 13,3 Pa (O.J. Tore), The aamc 
flexlMe cufihlou formed by opeu-coU foam purRoU with nlfroRon in u«od to 
Hupport the KevXar aircraft Hkin. TIiIh liunilatlon nyHiem la 27.9 mm (1,1 In.) 
thinner than the closcd-'ccn foam deolRu and atXll ruHultM in Tcuh fuel bolloff. 
However, it would be a more exponalve ayatom. Kxperlmoutal teatluR and de- 
velopment ol both ayatema hua been recommended to expl ore the potent la J aUvau- 
tURoa and ahorteomlnga of each. 

Thu volumea forward and aft of each tank in the aircraft (aee 1 ij’urea 1 
and 2) arc intended to be actively fluahed with air whenever hydroRen (or 
methane) acnaora located in the apace detect leaked fuel. In flight tliu air 
purge can be provided by opening vonta and uaing ram air. On the ground the 
poaitlve air flow can be provided by blowera inatulled for the purpoae. 

A unique potential problem poaed by the location of the tanka for the 
cryogenic fuels in the fuselage la the question that surviving passengers and 
crew could be faced with the hazard of expoaure to the fuel itself. Very 
brief contact with LH^ produces the effect of a burn on human skin. Exposure 
of large areas of the body to cryogenic temperatures for even brief periods, 
of course, would be fatal. 

It Is always possible to postulate circumstances Ju which large quantities 
of fuel are spilled and some of it drains into the passenger compartment. 
However, in view of the design of the tank and Its installation it is con- 
sidered highly improbable that there would be a survivable crash where this 
would happen. To illustrate, two pressure bulkheads separate the fuel from 
the passengers. As previously pointed out, the cryogenic fuel tanks are de- 
signed as pressure vessels to an operating pressure of (21 psia). Another 
pressure bulkhead forms the ends of the passenger compartment adjacent to 
both the forward and the aft fuel tanks. Both the end of the fuel tank and 
the adjacent bulkhead leading to the passenger compartment would have tu be 
ruptured to permit direct access of fuel to the passengers. Tu explore this 
eventuality an experimental program has been suggested (see Section 7, Recom- 
mendations for Future Work) wliich would involve crashes of suitably modified 
surplus aircraft. 
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3.3.2 raooutu'm: cublu wa.U auructuiu. '• One ol thts ImiKutum. t;auk'i wf thi.« 
study was (o uvaluatu thu rulativo ha>:m.'d ol paHuuuB'si'fJ Inoidc the oatjlu wltuu 
thu outaidu ol tlio luaolu^.p la OKpooiul to flj’ti from Hi)mod fuoi. To 
ovttluutu this lm;5ftrd it lo uottoouary to know the Htniotura] dooiau ol u 
tyjJitiaJ. Mcotion ol the luutslago mid the mutei'ialM involved in its eoMsUuctloii 
so the thorituil iiroiuirtioo emi he determined. 

View of ligure H iu u repruseiitutiva pieture of how the sidewall of 
a moderu iKisseiigur trmuijyort is constructed. The Lockheed L-lOU design is 
used us liluot ration. In the section where the windows are located the 
aircraft skin is reinforced with a doubler. J,'’ruraoa are located 20 inches on 
center. Above and below this window belt is a section where fail-safe straps 
of titanium are employed to reinforce the aluminum skin. A cross section 
through the wall in this area of the fuselage is illustrated j,n figure 9 with 
material callouts noted. 

The aircraft akin is clad 2024-T3 aluminum of 2.29 mm (0.09 inch) thick- 
ness. A layer of AA fiberglass insulation 57.2 mm (2.25 in.) thick is retained 
about 12.7 mm (0.5 in.) from the inside surface of the skin. In turn, spaced 
about 9.5 mm (3/8 in.) from the inside surface of the fiberglass insulation is 
a honeycomb panel which serves as the Interior wall of the cabin. The honey- 
comb panel consists of a core of crushed iioraex with surface sheets of hard 
fiberglass on both sides. The total thickness of the honeycomb panel is 
3.18 mm (0.125 in.). 

A sectiun showing the construction of a window Is presented in figure 10. 
Ihe 825 cm (128 in ) clear area consists of three panels, the outer two 
being structural members, the innermost being simply a protective layer to 
prevent scratching of the more critical panes. It also serves a function of 
permitting circulation of warm dry air to prevent fogging. 

The thermal characteristics of the cabin wall section are presented in 
section 5.5, along with an analysis of heat transfer through the wall as a 
function of time as a result of exposure to various fire conditions. 
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I Jh 






L-1011 cabin window aanurably. 


3.4 Aircraft Accident Records 

The sources of accident data used in this review were; 

a. National Transportation Safety Board (NTSB) 

b. luternationnl Civil Aviation Organiaation (ICAO) 

c. World Airline Accident Summary -- Civil Avi . inn Authority (CAA) 

d. Lockheed Company Accident Riles 

e. Miscellaneous Reports, Articles. 

The NTSB records, particularly accident suiimiari.ef:i (references 15 and 10) 
for the years 1964-1977 , were the primary source of information for compiling 
statistics. The other data sources, us well as mure curt' nt accident reports, 
were used to augjnent the NTSB records. A complete evaluation of the accident 
data is provided in reference 17. 
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Tabic 3 lii’ovldcs a auimnavy oi 783 NTSB accidcnta for tho porlod iy64--iy77 
uu fuiuaiouH of injury index and oiHu-auionai mode. Thin total wan reduced to 
341 since numy of the accidento were not pertinent for the purponen of tliio 
study. Those accidents eliminated from further contiideration were: 

• Small aircraft (TOGW < 5670 kg (12,500 lb)) 




Nondesimi-relatcd items (turbulence, hail. 

Improper safety practices (unfastened seat 
aisle injuries) 


sabotage) 

belt, passengur/attendant 


• Unrealistic Impact conditions (in~air aircraft collisons, collision 
with on-ground personnel) 

• Miscellaneous accidents (inadvertent door opening, static operations). 


The reduced number of accidents is shown in table 4. They consist of 
accidents due to: 


• Controlled and uncontrolled collisions with ground/water 

• Undershoot 

• Stall 

• Hard landing 

• Wheels up, retracted gear, gear collapse 

• Overshoot 

• Swerve 

• Obstacle impact or colir'sion. 

Figure 11 shows a comparison of the number of occurrences of fatal, 
serious injury and minor/non-injurious accidents for the original 783 data 
set and the reduced 341 accident data set. As can be noted in figure 11, 
both sets of data are consistent with each other in that both show uncontrolled 
and controlled collisions with ground resulting In .? substantial number of 
fatal and serious injury accidents; while stalls, ovo.r^dioots, uaderahoots and 
accidents involving obstacles, result lu a moderate number of fatal and serious 
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iAJiU'. i. ijUMMAUY Ol-' ACCIDENT TYl’Kti JiY QPliUA'nONAI. MODli 


t9C4 77 N rSB Oato 


Primary Acqidnnt Type 


I. Severn Impact 

A. Cumralled collisiop 2 0 Q Is 4 1 

B. Uncontrolled Collision 2 117 10 

C. Undershoot 4 2 20 

P Stall 1 2 8 2 1 0 


Taxi Flight S 
XY ZIXYZ XY 2X 


II. Moderate-High Sink Speed 

A. Hard landing 

B. Gear Collapse 

C. Wheels up 

0. Retracted Gear 

E. Swerve 

F. Overshoot 

III. System Malfunction 
A. Engine Malfunction 


2 1 27 

1 1 4 0 2 25 

0 0 10 1 18 

0 0 3 0 0 11 

1 1 7 0 0 25 

4 6 16 


1 9 11 0 4 3 


B. Prop/rotor malfunction 2 0 2 1 0 1 

C. Airframe failure 0 5 1 0 0 3 

0. Fire/explosion 0 2 2 0 1 4 


IV. Collision With 
A. Trees 


5 0 3 


B. Ditches, fence, seawall 1 3 11 0 0 


C. App. lights, wires 


0 1 2 0 2 4 


D. Obstaclos (bldg., auto) 10 3 1 0 9 

E. People 2 0 0 1 1 0 

F. Aircraft 2 0 0 0 1 0 


V. Miscellaneous unknowns 0 2 5 1 4 8 

VI. Turbulence and misc. flight 10 10 4 1 

occurrences 


9 1 0 

11 0 0 

1 0 0 


1 2 0 


2 3 20 4 6 0 

3 0 1 4 6 0 

0 3 4 0 5 

33 3 0 70 55 


1 0 1 

0 0 1 

0 9 0 0 2 0 0 1 

0 0 too 

0 17 10 1 10 0 2 

41 3 15 2 1 18 2 1 1 

1 203 6 


17 27 C? 44 35 179 4 10 46 48 223 46 7 31 10 1 1 0 

96 258 62 317 48 2 


X « No. of accidents Involving fatalities 
Y ~ No. of accidents In which highest injury index is severe injury 

Z ==< Nu. of accidents In which o .y minur/nu injuries ware sustained 

(T) « Total no. of accidents 


(T| 

Totals 
X Y 

‘1 

4 . 

196) 

54 


30 

(32) 

26 

5 

1 

(23) 

20 

2 

1 

(26) 

4 

2 

20 

(16) 

4 

3 

8 

(176) 

9 

ii 

152 

(30) 

2 

1 

27 

(47) 

2 

6 

40 

(20) 

0 

1 

19 

(16) 

0 

0 

16 

(37) 

1 

2 

34 

(26) 

4 

6 

16 

(131) 

20 

il 

171 

(54) 

3 

17 

34 

(21) 

10 

7 

4 

(21) 

4 

5 

12 

(35) 

3 

11 

21 

(95) 

30 

_9_ 

66 

(10) 

6 

0 

4 

(6) 

2 

3 

1 

(10) 

0 

3 

7 

(26) 

2 

0 

24 

(6) 

5 

1 

0 

(37) 

16 

2 

20 

(66) 

6 

44 

18 

(217) 

2 207 

8 


121 327 335 
(783) 






















TAULK A. 


UJvl)UCJ-:i) SUMMARY (U’ ACCIURNT TYPMS HY OPUUA'i.'I,ONAj, MOUK 
UASJ51) ON KiaJUOJi)) NUMUJ-;H or ACCIDRNT CANUIDATI.'IS 


1964-77 NTSB Uate 


Priniarv Acclilont Typo 

I. Severe Impact 

A. Controlled collision 

B. Uncontrolled collision 

C. Undershoot 

D. Stall 

II. Moderate-High Sink Speed 

A. Hard landing 

B. Gear Collapse 

C. Wheels up 

0. Retracted gear 
{:. Swerve 
F. Overshoot 

III. System Malfunction 

A. Engine malfunction 

B. Prop/rotor malfunction 

C. Airframe failure 
0. Fire/explosion 

IV. Collision With 

A. Trees 

B. Ditches, fence, seawall 

C. App. lights, wires 

D. Obstacles (bldg., auto) 


Takeoff 
X Y Z 


0 0 16 

1 1 2 

4 

1 4 2 


Landing 
Y Z 


4 1 

1 0 

2 IB 

0 0 


1 24 

2 23 0 0 

0 IS 

0 10 0 0 

0 16 0 1 

6 1G 


2 2 
0 1 

0 3 0 0 

14 0 1 


0 3 
0 0 
2 4 

0 2 0 0 


3 1 0 7 (191 


0 0 
6 0 0 


X = No. of accidents involving fatalities 
Y = No. of accidents in which highest injury Index is severe injury 
Z " Nu. of accidents in which only minor/no injuries were sustained 


1 ( 10 ) 
1 ( 12 ) 


11 18 35 37 21 144 0 2 26 15 1 31 

64 202 28 47 


37 9 


(28) 

23 

4 

1 

(11) 

8 

2 

1 

(14) 

4 

2 

18 

(7) 

2 

1 

4 

154) 

8 

12 

134 

(27) 

2 

1 

24 

(42) 

1 

3 

38 

(16) 

0 

0 

16 

(14) 

0 

0 

14 

(29) 

1 

2 

26 

(26) 

4 

6 

16 

(80) 

8 

15 

57 

(38) 

2 

9 

27 

(7) 

4 

n 

3 

(16) 

1 

3 

12 

(19) 

1 

3 

15 

(37) 

10 

6 

21 

(10) 

6 

0 

4 

(5) 

2 

3 

0 


63 42 236 

(341) 


(T) ~ Total Nu. of accident. 
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Number of accidents 
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Injury aculduntH un or mrountl HirportUt Couvoruoly, wheola-upi rotractod goav i 
sworvo, hard lauding and goar collapun ruuull; in low latal *>r aurloim injury 
acnldniita . 

Thn NXUJJ data rcpmannlH lowH than 2‘1 pornnut ol the ti>tai world anoldonta 
during the V)bh"Wn perijd. I'Mgnre 12 ahowii a cuiiiipar l non t)l’ the lo.dnc.od 
341 Ni'SB accidents and worldwide accident types. Xhu trend between both 
accident seta ol data indicates that the reduced NTBll accident data is 
representative ol the worldwide accident history Insolar as severity ol 
accident types are conccrued. Table 5 suminariises the comparison between the 
original set ol NTSU data, the reduced set ol NTSB data, and the worldwide 
data set. 

A total of 108 available accident reports were obtained from the reduced 
NTSB data set. As part of the study described iu reference 17 these 
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CC - CONTROLLED COLLISIOSW/GnOUNO 
UC - UNCONTHOLLEO COLLISION W/GROUND 
US - UNDERSHOOT 
HL - HARO LANDING 
WU - WHEELS UP 
RG - RETRACTED GEAR 
QC ~ GEAR COLLAPSE 
OS - OVERSHOOT 
OttST - OBSTACLE 
•THRU'mAHCH 1379 



^ 2707 ACCIDENTS ^ 

WUHLOWlOE 1964 79* 

Figure 12. - Comparison of NTSB and worldwide accident summary data. 


341 ACCIDENTS 
NTSB 1964 77 
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TAJlLli I). 


COMl'ARiSON OR RATAL ACCimiNT RliUCKNTAOlSS ROK 
NTUU ANU WORLDWIDE ACCIDENT EUMMARIES 



NTSB1964-77 
(7B3 Accidents! 

NT8B 2964 77 
(341 Accidents! 

Worldwide (1964 - 79!* 
(2707 Accidents! 

Accident Tyite 

No. Fetal 
Accidents 

No. Total 
Accidents 

Fatal: 

No. I'atal 
Accidents 

No. Total 
Accidents 

% Fatal: 

No. Fatal 
Accidents 

No, Total 
Accidents 

% Fatal: 

Cantrolled collitlon 

26 

32 

81.3 

23 

28 

82.1 

58 

100 

58 

Uncontrolled collision 

20 

23 

87 

8 

11 

72.7 

27 

40 

67.5 

Stall 

4 

15 

26.7 

2 

7 

28.6 

16 

33 

48.5 

Undershoot 

4 

26 

16.4 

4 

24 

16.7 

14 

37 

37.8 

Hard landing 

2 

30 

6.7 

2 

27 

7.4 

2 

51 

3.9 

Wheels up 

0 

20 

0 

0 

16 

0 

0 

27 

0 

Retracted gear 

0 

16 

0 

0 

14 

0 

0 

67 

0 

Gear collapse 

2 

47 

4.3 

1 

42 

2.4 

4 

152 

2.6 

Swerve 

1 

37 

2.7 

1 

29 

3.4 

8 

88 

9 

Overshoot 

4 

46 

8.7 

4 

26 

16.4 

4 

76 

5.3 

Collision with obstacle 

10 

62 

19.2 

10 

37 

27 

16 

51 

31.4 

•Thru March 1973 


accideuCs ware clasaiflad in accordance with the accident types noted in 
figures 10 and 11. The accident types were then classified into three crash 
situation categories. These crash situations are defined below asj 

CCOi Cround-to~Cround Overrun Condition (table 6) 

AGUL: Air-to-Cround Hard Landing (table 7) 

ACl: Alr-to-Ground Impact (table d) 

ihirty-two of the 108 accidents were considered as "impact nonsurvivable" 
or "wipe-outs"; that is, everyone onboard perished. Since Scenarios 2 and 3 
in the subject study are concerned only with crash fire hazard in otherwise 
survlvable accidents those 32 accidents were eliminated from further consider- 
ation. In addition, ten miscellaneous accidents, which yielded little or no 
useful information as to causal data were eliminated from the accident files 
leaving a total of 66 accidents available for data analysis in the three 
selected crash situations: GGO, AGHL, AGt. 
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TABLE 6. - GROUND-TO-GROUKD, OVERRUN CONDITION, GGO 



TABLE 7. - AIR-TO-GROUND, HARD LANDING, AGHL 


I 
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dij U’).fia dui iiut iui,); uJ. whiil. tumwLiuattti n aurvivati-Lo acciUaiu:. Tli*.' 
N.iUiuual Ti'iuuiinii’t.ai.J.ou tjataiy liimrU (NTl>iJl) duU.atifj a aiirvivaliXu ai.a'.idt!iu. 


a« oue iu vjhiah tlui oaaupiably vulumo iu cmiMidui’ud Hul'liainiiilly Intact folluw" 
iny a ctaali to have ptovidud ptotootion tot tla' oucupautu and tho loads aru 
doomed to have Leou within tolurabXu limlta. Thus, by NTSU rules, aocldonts 
which some pooplo aurvlvu are somutimus aw nou'-survivablo and cou" 

vursuly, accidonta which produce no survivors arc Hometlmca labeled survlvablu. 
On the other hand, the I’cuural Aviation Authority (FAA) employs a straight- 
forward dulinition which rtaults in a classification of survivable if at least 


one person survives a crush. For the present study wu arc using the FAA 
definition. 


Again as part of the study described in reference 17, the accident data 
was then assembled in accordance with the type of airplane engine arrangement 
eraployeu! wing pylon-mounted engines, aft fuselage-mounted engines, and 
combined alt f uselage/wing pylon-mounted engines. The distribution of accidents 
as they relate to these different engine mounting arrangements, is shown in 
table y. Of the total of 66 accidents, 30 hud an engine mounted in or on the 
rear of the fuselage and thus had fuel lines running axially through the 
fuselage which were subject to damage from fuselage breakage/separation. 


TABLE 9. - DISTRIBUTION OF ACCIDENT DATA AS RELATED 
I'O AIRPLANE ENGINE ARISING EMENT 



Candidete Crash Situation* 



GGO 

AGHL 

1 AG1 


Engine 

Arrangement 

GG1 







Landing 

Takeoff 

GG2 

AG1 

AQ2 

AG3 

AG4 

Total* 

A** r 'iisge 

4 

2 

4 

7 

2 

3 

3 

25 

t. 

6 

6 

6 

3 

1 

4 

11 

36 

One Att Fuulage - 
Two Wing Pylon 

- 

2 

- 

1 

1 

1 


6 

Totilt 

10 

10 

8 

11 

4 

8 

14 

66 

1 *0ifln«itlnTibl«6, 7, tnd8 





laUitui IQ, il auU 12 h1m>w l-tif Ci u(|i«!iicy i»l uanun’Muet* ut fm.*! apl J lajjii lit 
ittlauiuu Uu Lilt! <icuu»'it!iu;t! <jl luiittlajn! HujiaraLloa fur l1u> tIillurauL uuniiu! 
rtmUi^urtiU aijplautia. Aa ran ba aultid by Lbu data in tliuau Labiuti iutarlaaM 
iiu|»a) al-itiu/bi’i'ak urc.urfi al wr I I drlinrd liai'd iioJut.H, ,l.Ut, winy, Itiadiufs 
r'dyr winy, Lr.iiliiiy, t!dy,i,' (T.li.), Inrward nl niiiiH'.nnay.t! , and aJ I, uj 


rab in . 
only nun 


in Lilt' tU'r J.dnn Ln involving Llin 30 airnrait. wiLb ai L'^tiKHinLud nn^ylnna 
iniitaiuai o) iuanlai'.t! aeparat inn wan unuiid an a iiunaib.le (.mu, at' nl 


I iiu.l. api I I , 


Tabica 13 and lA nlinw Llin rulutivu iirn hazaid and 


atructural ayatum 


danaiii' invnlvamonL Jur al l tbu aiicral t aa a iunctioii oi their 
mei;t, Ul: a Lntal nl 3b .leeidenta involvlii}' alrplanea with win|< 


engine arrange- 
pylon-mounted 


unginen, 2^ acxidents allow lire occur renee.s or the potential of a lire hnaard 


truii fuel apillage. Wing fuel tank rupture occurred an the stated cause in 10 
ot the 24 lire hazard accidents, lollowed by 8 oeeurrences ol wing iuel line 
aeve/ence. Misceliuneuus or unstated causes account lor 6 lire hazard 
ac-cid Jilts. Generally the iiiiscelluneuus cause is associated with a wing root 
separation or wing tank rupture. 


TAliLii 10. 


lUiLATiUNSUil* liETWKliN 


FUSELAGE BREAK AND FUEL SPILL OCCURRENCE 


- GGO GK. 11 SITUATION - 


Engine 

Arrangement 

Number of 
Accidents 
Involved 

Occurrence of Fuselage 
Break and Fuel Spill 

Occurrence of Fuel 
Spill end No Fuselage 
Break 

Occurrence of Fuselage 
Break and No Fuel Spill 

Aft fuselage 

7 


1 

,.5 

Wing/pylon 

9 


4 

2 ’ 1 ' 

One aft fuselage/ 
two wing pylon 

2 


2 



(^13 Break uccurreil at wing L.E. 

( 2 ) Breaks occurred at wing L.E. and T.E. 

(^3^ Break occurred furwar:' of emponnage, fuel spill occurred from wing fuel tank. 

( 4 ) Breaks occurred at wing L.E. and T.E., fuel spill occurred at wing root. 

(53 Break occurred forward of empennage. 

(§) Breaks occurred at wing L.E. and T.E., fuel spill from wing fuel tank or engine fuel line. 

Breaks occurred at wing L.E. and T.E., fuel spill frum wing fuel tank. 

(B) Break occurred at wing L.E., fuel spill from wing fuel lank. 

Break occurred at wing L.E., unknown leak source. 
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i'AUJ.li il. ~ KhLATiUNtiim* liK'i'WliliN FUaiiLAGIi UKlvAK ANU l-'UI-l, Sl'XLL OCdURUliNCfc: 

- AUIIL (JHASII aXXUATXON - 


^ Break occurred forward of empennage. 

Breaks occurred aft of cabin and forward of empennage. 

^ Break occurred at wing T.E. 

Break occurred forward of empennage, full spill occurred at wing root. 


Engine 

ArrangHinent 

Number of 
Accidents 
Involved 

Occurrence of Euselago 
Break and Fuel Spill 

Occurrence of Fuel 
Spill and Nq Fuselage 
Break 

Occurrence of Fuselage 
Break and No Fuel Spill 

Aft fuselage 

S 

1 ® 

1 

1 0 J 1 ( 2 ) ^ ( 3 ) 

Wtng/pylon 

2 


1 

1© 

One aft fuselage/ 
two wing pylon 

1 


1 

- 


rAbLE 12. - RELATIONSHIP BETWEEN FUSELAGE BREAK AND FUEL SPILL OCCUiytENCE 

-• AGI CRASH SITUATION - 


Engine 

Arrangement 

Nomber of 
Accidents 
Involved 

Occurrence of Fuselage 
Break and Fuel Spill 

Aft fuselage 

8 

1©4©1© 

Wing/pylon 

18 

101012® 

One aft fuselage/ 
two wing pylon 

2 

© 

© 


Occurrence of Puel 
Spill and No Fuselage 
Break 


Occurrence of Fuselage 
Break and No Fuel Spill 



(p Breaks occurred aft of cabin and at wing L.E, and T.E., fuel spill due to engine line or fuselage separation. 

Fuselage destroyed or broke up into several sections and fuel spill attributed to wing tank rupture. 

@ Breakup occurred forward of empennoge, fuel spill unknown, could be fuel tank or fuel line rupture, 

0 Breakup occurred forward of wing L.E., fuel spill from winy tank. 

V Breakup occurred forward of empennage, fuel spill due to engine separation. 

0 Breakup occurred at wing T.E., fuel spill due to engine separation. 

© Aircraft were destroyed in 4 accidents, breakup occurred aft of cabin 3 times, forward of empennage 2 times 
aft of cabin and forward of empennage 2 times and at wing L.E. once. All spills are associated with fuel tank’ 
rupture, engi ne separation, or wing separation except in 2 cases where spillage cause is unknown 










TAUUi J3, ■ I’lHl'; UAl^.ARU tJUMMAKY, A1<J, AJUl'J.ANI'; UONJ’JCiUlUVTlONii 


AU1 I A02 I A«a 1 A(H I Tumi 


Nuinlier of Accidonti 

FIHUHAZARI] 

Total Number Firut Occurred 

Fuel Cuuteinmont, Struc. Related 
(Fuel SdIII) 

Friction, Ground, Hydraulic Line 
Related 

Unstated 

Fuel Spill, No Fire 

Total Number of Fires and 
Potential Fires 

ESTIMATED CONTRIBUTIONS TO 
EIRE HAZARDSfn 


• Wing Fuel Line Severed 

a. Englne/Pylon Separation 

b. Wing Root/Inboard Failure 

c. Wing Outboard Failure 

d. L.G. Penetration 

e. WlngL.E, 

f. Miscellaneous 

a Fuselage Fuel Line Severed 

a. Fuselage Soparatlen 

b. Engine Separation 

c. Landing Gear Penetration 
■ Wing Fuel Tank Rupture 

a. L.G. Collapie/Penotration 

b. Columnar Impact CO 

c. Contour Impact® 

d. Frontal Obstruction ® 

e. L.G. Collapse and Fus. Impact 

f. L.G. Trunnion Attach Fitting 

g. Miscellaneous 


* Refined in Tables G, 7, end 8 
® Pole, Tree, Pier, Light 
@ Ditch, Ravine, Hilt, Embankment 
@ Fence, Building, Wall 

^ Incluoas possibility of either or both occurrences on 
same accident 

Fuselage Llne,/routing designed as a result of accident 
in 1966 


I Wes') 


(*X9) 


^ Obstacles Involved, svjng separation occurred 
m Fuel cuntainment related fires and spills 
@ Nose gear collapse, severed hydraulic line 
Possible fuel tank rupture 
(10) Wing Root Area 


1 


TAllLK 14. - aTRUClTUltAL DAWAGK SUMMAUV, AU. AXRPLANJi CONPiaUliATIONS 



Candidate Crash Situations* 


G(I1 

flG2 

aq; 

AG2 

AQ3 

AG4 

Total 

Number u( Aecldenti 

20 

B 

11 

4 

a 

14 

GO 

('uieliife Bieak/Sennretlun 








e. Aft of Ceukplt Fwd, ruiolege 






3 

3 

b, Winy, Lending tldgo (L.ti.) 

2 

1 




1 

4 

c. Winy, Trnlling Edge (T.li.) 



2 

1 

1 

1 

4 

d, forward gf Bmpeiinayo 

1 

1 

2 

3 

e 

G 

tt. Wing, LE. end T.Ei. 
f Aft of luck|jlt and f'wd, 

2 

3 

1 


2 


2 

of Empennage 







g. Thren . • T'lur Breaks 


2 



2 

1 

4 

G 

h. r^ selii';.? flutroyed 




1 

2 

7 

Engine Separation 








Fuselage Mounted Configuration 

2 

2 

3 

2 

1 

4 

2 

12 

Wbig Mounted Configuration 

6 

1 

•s 

L 

1 

S 

20 

Fuselage and Winy Mounted 

2 


1 

1 

1 

6 

Wing Failures 








Root-Inboard 

2 

3 


2 

a 

11 

1 

2 

21 

9 

6 

Outboard-Tip 

2 

3 


2 

1 

1 

Mitcellaneous/Unknawn 

2 

1 



landing Gear Collapse 








MIG Collapse Only 

4 

2 

6 



3 

15 

7 

NLG Collapse Only 

5 

1 

1 




MLG and NLG Collapse 

7 

3 

1 

3 

3 

3 

20 

Occurrences of; 








Floor Buckling 

1 

3 

4 

3 

1 

1 

3 

15 
8 

16 
11 
13 

Cockpit Crushing 

1 

1 



e 

Lower Fuselage Crushing/Abrasion 

5 

3 

4 

2 

9 

Overhead Racks/Panel F Allures 

3 

4 

3 



1 

Q 

Jammed and/or Deformed Doors 

b 

3 

1 

1 


Slide Related Failure 

3 

1 



Galley Equipment Failures 


1 




£ 

b 

1 

6 

Tnc/Wheel Oamago 

e 

1 





‘Defined in Tables b, 7. and B 





' 



Numbers in table indicate number of accidents in which stated events occurred 







Ui a uotai ol' at-'i’-iduntu involving ni u- monnteU engine aivi»laneH> 

.1!) iieeitleniy nhow lite oeevuienee , n: ilre liotential firoiii a iuel spill , in 
Uliene 11) lire haj'uj.’cl aecldiuitn , aii iutieJage lu’eak up oeeurreU 12 iiiiiea. 
lU)weve).’, in only one ol iliene 12 -.nstaneea la it nuteU Utau the iuaelagi! 
scparatlon/hreak may have led to the i'uel spill. Again wing fuel tank ruptuve 
ia a prevalent occurrence (8)j followed by fuaelage fuel line aevorencu, 
usually engine separation (4)i and wing fuel line aeverencui particularly at 
the ruot (3). 

The explanation of the cauae of fuel leakage or fire ia oft timea 
unclear in the accident reportn. I'or example, in table 13 it is reported that 
in seven Instunees the eause of fire was unstated in the NTSli reports. Since 
speeifie causes of fires wliieh were not related to fuel spills arc identified 
(three eases) it can be presumed that the sevei’ unstated eases involve fuel 
spills but that the source uf the spill was in question. On the basis of this 
conjecture, it would be valid to conclude that 82 percent uf the time fuel 
is spilled, fire results. On the other hand, if the sever unstated causes uf 
fires are divided in the same ratio as the number of fires known to be due 
to fuel spills (31) to fires known to be due to nonfuel causes (3), then 
only 6.3 of the seven unstated fires would be due to fuel spills and the 
above percentage would drop to 81.2 percent. It therefore seems appropriate 
to eouelude that in more than 80 percent of the time when fuel is spilled 
in an aircraft crash, fire results when the fuel is Jet A. Further, since 
all the other fuels being considered in this study will ignite at least as 
readily as Jet A, the statement could be expanded by saying it would be true 
regardless of the fuel. 

Table IS summarizes the fuel containment related fire potential and 
structural system damage data. Fuselage break, engine separation, wing 
failure, and landing gear collapse occurrences, in terms of relative frequency, 
are similar fur wing pylon or aft fuselage engine-mounted configurations. 

While some percentages differ for the two different configurations, they a.‘e 
generally within ±10 percent of the total of all 66 accidents. 
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i'AlSl.K 


tiLlMl-UKV Ul' CONlAiNMIiNT UAJ'.AHD.S ANU liTKlUri'IJlt^M. bA.'UM.lv 

Aa A ^'u^'c:TI0 N OF aruanofmknt (;o^;l•'J.;ukA■^o^i 



Lngina ArrangiimQiit | 


Wing Pylon 

Aft Fusulogo 

All 

Number of Accidents 

36 

26 

66 

Fire Potential Fuel Containment Related 

18 (601® 

16(60)® 


Wing Fuel Lino Severance ® 

8 (22.2) 

Sd.T) 


Fuselage Fuel Lino Severence 

- 

4(16)® 


Wing Fuel Tank Rupture ® 

10(27.81 

8 (32) 


Structural System Damage CD 

1 



Fuselage Break/Separation 

20 (55.6) 

16 (64) 

39 (69.1) 

Engine Separation 

20 (55.6) 

12 (48) 

37 (66.0) 

Wing Failures 

20 (55.6) 

11 (44) 

36 (54.5) 

Root-Inboard 

13(3G.1) 

5 (20) 

21 (31.8) 

Outboard Tip 

5(13.9) 

4(16) 

9(13.6) 

Miscellaneous 

2 (5.6) 

2 (8) 

6(9.1) 

Landing Gear Collapse 

22 (61.1) 

16 (BO) 

42 (13.6) 


( ) Denotes Fuel Conteinment and Structural Damage Data hi percentages: based on total number of accidents for 
the particular engine arrangement configurations noted at the top of the respective columns. 

(T) More than one system can be damaged in an accident, therefore, percentages can total more than 1Q0%. 

(2^ Seven fire potentials noted as wing fuel line severence could also be wing fuel tank ruptures. 

( 3 ) Only one due to fuselage separation. 

( 4 ) Does not include 1 fire occurren''e in which containment relationship is unclear. 

(53 Does not include 6 fire occurrences in which containment relationship is unclear. 



Thn location of the enRlnes on the wlr.R, or on the aft fuficlaRe. does not 

algnlf leant ly alter the accident statistics with rcRard to overall structure 
participation. 


A detailed evaluation of the crash conditions postulated for the suh.lect 
study Indicates that the three crash situations and their subsets (;c2, 

AGl, A(t 2, AG3, AG^), are applicable to crash scenarios 2 and 3 as follows! 

Crash scenario No. 2, defined as follows! "A survivablo crush lundlnc, 
or failed take-off. where damage to fuel tankigc or lines results in massive 
release of liquid hydrogen after the aircraft has come to rest," would result 
in damage similar to crash sitaations GGl and GG2, which are essentially 
ground-to-ground accidents occurring on the runway or as a result of over- 
running the runway. The accidents in these scenarios involve airplanes running 
erf the runway at a speed in the vicinity of 80 kts and impacting terrain, 
or obstacles, at lower speeds; e.g., 40 kts. On an average, the airplane 
comes to rest within 300 meters past the runway threshold. For these types 
of accidents, 14 of the 29 accidents involved fuel leakage. Wing tank 
ruptures account for 11 of the 14 spills distributed as follows; 

7 contour, columnar, or frontal impacts* 

4 landing gear penetration 

Two of the 14 spills were associated with wing engine/pylon fuel line breaks 
and one of the spills was associated with a fuselage fuel line break (could 
be either a fuselage break or engine separation). 


Fuel spill from only one side of the aircraft is characteristic of this 

type of accident. Fuselage breaks are generally not severe and the separation 
usually isn’t complete. 


Crash scenario 3, 
failed take-off, where 


defined as follows; "A survivable crash landing, or 
damage to fuel tankage or lines results in massive 


release of liquid hydrogen upon impact and during aircraft deceleration", 
would produce damage similar to AGl, AG2, AG3 and AG4 crash situations. The 


*3 of these spills could be wing root /inboard failures. 
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,U):"i-()''tjrtnmd in uIh’.bu iiivolvt'. high inrward and sink 

spuads. A iorward spaad of at .l.oast 120 kt s was ausumad for tha purpose, of 
uauiinauiug aradikla dainaga .laadiug to fuel ap . ilaga rates and quantities. 

AOl and A(.i2 aacidents oeaur in the proximity of the airport, while A03 
and AGA ae.eidents oft-times oaeur substantial distances from airports and in 
hostile surroundinas. For these types of crashes, 24 of the 37 a.iidents 
involved fuel leakaac. Wing fuel tank ruptures accounted for 18 spills. In 
the majority of the accidents, columnar (tree), contour (embankment, hill) 
and frontal (house) impacts caused failure of the wing either inboard at the 
wing rout or outboard of the nacelle. Fuel line leakage occurred In two 
instances due to wing engine pylon separation, two Instances of fuselage 
engine pylon separation, and two instances involving rupture of a line in 
the wing leading edge. 

Fuel spillage from both sides of the aircraft Is characteristic of this 
type of accident due to the density of obstacles in the crash environment. 

The severity of the crash conditions for this situation is such that complete 
separation of fuselage sections on break-up can be expected. 

In summary, in the accidents studied It was found Chat although fuselage 
breakage/separation occurred just as frequently as did wing damage, the fuse- 
lage damage did not usually result in fuel spillage, whereas wing damage did. 
Further, in the case of cryogenic fuels, it must be recognized that they will 
be contained in tanks designed to an ultimate pressure significantly greater 
(30.2 psi vs. 18 psi) than that of the rest of the fuselage In which they are 
mounted. Therefore, it can be concluded that ';be fuselage-mounted tanks of 
the cryogenic fuels will be more likely to survive a crash than the rest of the 
fuselage, i.e. , breakage will occur In the weaker sections. This should lead 
to markedly better survival statistics for the cryogenic fuels. Fuel spillage 
will occur, generally, only if fuselage sections separate so far that lines 
are pulled apart. Engine separation from the aircraft, whether wing or 
fuselage -mounted, will be the principal source of fuel spillage for the 
cryogenic-fueled aircraft. In contrast, the conventionally fueled aircraft 
will continue to have both this problem, plus the vulnerability of wing tanks 
which has been shown, in the foregoing analysis, to be the major factor leading 
to fuel spillage. 
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The ocher Cwo crash acennrlos, Nos, 1 and 4| rcprcHonced casca In which 
crash datnoRc wna noc rcIcvanC. 

Ccmih scenario No. ] , sCnCed as f'onows: "A normal landing;, or ground 
accident, which results In fuel system Insulatltjn damoRc and/or fuel system 
damnRO pcrmlttlnR liquid liydroRcn to vent, escape, leak, or run out of a 
punctured tank or broken line," should be considered an accident In which only 
small leakaRc rates could occur because control of flow via shutoff valves 
and boost pump ussRe are controllable by the crew. Tills situation would come 
about as a result of a lonR developlnR problem (metal fatigue) or. It could be 
maintenance related as opposed to crash damage. 

Crash scenario No. 4 Is a situation where the aircraft plunges to the 
ground and all on board are killed by the impact. The question to be resolved 
is, what is the hazard to the surrounding community represented by the in- 
stantaneous release of the fuel aboard at the moment of Impact. 


i.5 Scenario Damage Definition 

In light of the damage which is apt to be inflicted on the subject air- 
craft in the postulated scenarios, as described in the preceding section, an 
analysis was made to determine what fuel leakage might be expected to result 
with each of the four fuels. 


I by way of review, the cryogenic fuels (LH 2 and LCH^) are stored in 

, insulated tanks located botii forward and aft of the passenger compartment. 

I i’he non-cryogcnic fuels (Jet A and JP-4) are stored in the wing and in 

► '^®''ter section taiiku. As listed in table 2, the subject aircraft carry the 

I following total fuel loads: 


LH2 

LCll, 

4 

Jet A 


25,600 kg ( 56,460 lb) 
69,040 kg (152,200 lb) 
84,780 kg (186,900 lb) 


JP-4 

All of the aircraft have four 


84,780 kg (186,900 lb) 
engines mounted on pylons below the wing. 


Tlie location of the fuel tauks and engines dictates the primary effect 
of impact damage on fuel leakage. Wlieii fuel tanks are located in the wing, 
impact damage to fuel lines with.in the wing is of little concern because 


t 

i 

I 


( 

i 
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of the oxeesfllvo fuel tank damoRe that muat occur bofore the fuel Unea are 
Impacted, In all other arena, such aa cuRlnea, cuRlno pyloua, and fuel Tinea 
external to fuel tanka, fuel line leakaRe la of concern. 

liach engine la aaaumed to be supplied fuel from its own individual fuel 
tank. I'wo fuel pumps were used in cadi noncryogenic fuel tank and thtee in 
each cryogenic tank. 

Although as has been pointed out, the fuselage tanks of the cryogenic 
fuels are considered less apt to be seriously damaged so us to leak fuel in an 
impact survivable crash, compared with the wing tanks of the conventionally 
fueled aircraft, the problem at this juncture was not to assess probabilities 
but to decide on credible damage which could occur in order to aualyxe the 
consequences of fuel spillage. Accordingly, fuel tank damage was assumed to 
vary from a small puncture or a cracked weldment permitting a leak of 0.5 kg 
per second, to a major tank rupture. In the latter case, an average spill 
rate was determined by assuming that all of the fuel tanks in the designated 
area released their fuel within a 30-second period. 

Fuel line leakage is significantly affected by the operation of boost 
pumps. If the pump power source is turned off prior to the crash, spillage 
would be the result of gravity drain through a severed fuel line. For use in 
the subsequent analyses, minimum leakage was established by assuming one 
broken line only, fed by gravity or siphon action. Maximum leakage assumed 
all fuel lines in the area were broken with all boost pumps operating wlien 
applicable. 

Table 16 provides a summary of leak rates and quantities whicli were 
adopted to represent the LH^-fucled aircraft in the analyses of Scenarios 1, 

2, 3, and 4 described in the following sections. The fuel quantities listed 
in the table for Scenario 4 represent the amounts of hydrogen assumed to 
remain in the aircraft after, in the first case, a crash following taxi, 
take-off, and climb to about 1.8 kra (6000 ft). The second value listed is 
for the situation where the aircraft has flown its design mission and is 
near its destination when it crushes with only a little more tlian reserve 
fuel on board. 
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iAULt 16. " UltASU aiXUA'l’ION ANj) Sl'lLLACiJi l-’OR AlHCltAFl 


Situeilan 


Uainagn 


Scenario 1 Aircraft at rott. Smell intonial leak. Cracked wold on tank 


Scenario 2 Aircraft stopirod by Impact at 
':40kt$. Radial spill. 


Scenario 3 Damage occurs at 120 kts. Long 
axial spill as aircraft decelerates 


Engine loparation or fuselage lino break, 
a Boost Pumps Off 
• Boost Pumps On 
Tank pimctiired or ripped. 


LIU Spill 

Hutu Qe 
(kg/iec.| 


fSame as Scenario 2) 


Scenario 4 Cataitrophic crash. Large radial spill. All tanks ruptured. 



Corresponding leakage rates and spill quantities for the other fuels w« 
arbitrarily assumed to be in the ratios of the fuel quantities tanked, vi^, , 


69.040 
25,600 . 

84,780 

25,600 


2.7 for LCU^, and 


*= 3.3 for Jet A and JP-4. 


4. COm’UTEK MOUEI.B AND ^U^THODOT.O^.y 


The computer models and methodology used to calculate the spread » 
vaporisation, and dispersion of spilled aircraft fuels Irto the atmosphere 
are provided In this Sectloui For a fuel leak within the aircraft, a method 
Is provided for calculating the required air purge rate In the compartment 
to reduce the air/ fuel mixture ratio below its flammable limit. Finally, a 
thermal model Is described for calculating heat transfer rates and tempera- 
tures of an aircraft fuselage (both external and Internal) that Is exposed 
to flames from spilled liquid fuel. 

A.l Aircraft Fuel Liquid Spill, Spreading and Vaporization Models 

An analytical model was Initially developed to solve for the liquid 
spill, spreading and vaporization processes that would occur following the 
rupture of aircraft tanks or lines. Two .nodels were developed, one for axial 
spread of the liquid perpendicular to the direction of aircraft motion, and 
the other for radial spread of the liquid from a fixed aircraft position. A 
flat, horizontal spill surface was assumed, and no dikes or boundaries were 
present to limit the extent of the liquid spread. Concrete and soil (with 
10% moisture) were the two spill surfaces considered in these programs. 

These models are applicable to all liquids, including the cryogenic fuels 
and the conventional petroleum-based fuels considered in this study. 

The liquid spreading and vaporization processes were modelled using open 
channel, hydraulic flow methods '>f analysis. In this analysis, no combustion 
1s assumed to occur over the liquid spill and spread, and hence the unburned 
fuel vapors will disperse into the atmosphere. Applying both conservation of 
mass and horizontal momentum, two partial differential equations are obtained 
which define the liquid flow, as shown below for the axial spreading model 
case : 
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CQnscrvttClon of Masa: 


o|. pn,. *■: ■ 


— 1 




fK'W) 

fix 




0 


0 ) 


ConHcvvatlon of Momonuum: 


fit fix 




ili 

fix 


Y v|v|- 0 


( 2 ) 


whuro 

t • time from the start of the aplll 
X * axial distance from the center of the spill line 
V " mean horizontal velocity of the liquid flow 
6 ■ liquid flow depth 
g ■ acceleration of gravity 

Vg ■ rate of fall of still liquid surface* due to evaporation* etc. 
f " Fanning friction factor for open channel flow 

These partial differential equations were then converted to ordinary 
differential equations with respect to time* using the method of characteristics* 
as presented in Reference (18). These ordinary differential equations are 
then grouped In terms of the positive and negative characterisitcs of the 
^ slope, as shown below; 


Positive 

characteristic (3) 


Negative 

characteristic (4) 


. dV . . d6 ... 


vVg + Y VIVI- 0 


dx 

dt 


- V + \ 




dx 

dt 


- V ~ 


where \ " t yf%b ■ wave velocity. 
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Oh 


Computer HolutlnnH o£ tlip above differential ecjuatltms were then prof>;ranmiod 
for a Uiilvae UOO computer, ualnR mj explicit apeeified time interval method 
of solution, Blmllar to Section 15. A of Ueferonce (IB), 

in the above formulations, only the evaporation term Vj, and the. friction 
factor f are dependent on fluid propertioti. The evaporation term V., repre- 
aant.B the rate of fall, of a still liquid surface due to evaporation, bolHnp., 
or seepage of liquid into soil, etc. For cryoRonlc liquid boilinj^ on a sur- 
face, the evaporation t?rm can be found from the followinR rcl uttonship , 


V 


E 


P, h- 
L fg 


(5) 


where 

Pj^ “ liquid density 

h^g “ latent heat of vaporization 

q" ■ heat rate per unit area from the solid surface 

The transient heat flux q" from a semi-infinite spill surface can be predicted 
from Section 13-3 of Reference (19). For an initial time period, this heat 
rate is limited by the boiling burn-out heat flux of the spilled cryogenic 
liquid. After that, the heat flux depends on the solid surface thermal pro- 
perties, and varies Inversely with the square-root of time t. 

The term f represents a Fanning friction factor for open channel flow, 
derived from Section 4.2 of Reference (20), which must be expressed in the 
following form, 

*°'io [ifs * I 

^ n 

6 To 3 77 91 " parameter defining surface roughness. 

Manning's roughness coefficient, see References (20), (21), etc. 


where 

d - 

s 

n ■ 


58 


A6Vp, 

Ro " — " Reynolds number of the flow 
|j,j^ ■ liquid viscosity 

The first term in the logarltlun bracket represents the effect of surface 
roughnessi while the second represents the effects of flow Reynolds number 
for turbulent flows, The Ifannlng friction factor for laminar flow, f ■ 16/Re, 
was only used for low Reynolds number flows, where the laminar flow f exceeded 
the above turbulent flow f value. Hence, the flow friction factor and shear 
stress were defined for the full range of surface roughnesses and flow Reynolds 
numbers. 

In order to make a time Integration of the characteristic differencial 
equations, the boundary and Initial conditions must be specified. Figure 13 
shows the inlet flow ooundary and initial conditions at the spill location. 

Two inlet boundary conditions are possible, depending on whether the Froude 
number ratio of Inlet velocity Co wave velocity is less or greater than unity. 
At the Initial time t ■ 0, the spill surface Is completely dry, while the ini- 
tial inlet Froude number must be sl.O. 

The advancing free leading edge boundary condition is presented on Fig- 
ure 14. With Froude number Fr^l.O at the advancing leading edge of the flow, 
this fluid boundary cannot influence the method of characteristics solution 
interior to the liquid flow body. The total time derivatives of the depth 6 
and the flow per unit width (&V) were obtained from the partial differential 
equations, and expressed in terms of partial derivatives with respect to 
distance x. These equations both follow Che characteristic velocity ■■ V, 
which is the velocity of the advancing leading edge. 

A final boundary condition, that Involving receding free liquid edges in 
the final stages of fluid flow and vaporization, is shown on Figure 15. At a 
receding liquid edge, both the liquid depth 6 and velocity V are zero, and the 
momentum effects are small. The edge will recede due to vaporization, and the 
location of the edge can be computed by conservation of mass in the triangular 
wedge shaped edge elements. Receding edges can develop at a leading edge that 
stops advancing, or at the liquid Inlet position, x - 0, after the Inlet liquid 
flow has ceased. 
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I'roudo nurolMir TM ’ • 1,0 ut tlio lUlvuMuiitu lom1ln« mluit, 

Coinurving miiM iiiul inomontum, tim iiiutlul tiitforuiitiul o(|uutlmii nm uiml 


to form thii Ithullng »Hgt) modol: 

tl6 ()6 Oft 

— ■= ^ V -V„ 

dt <)t i)x 


i) I ft V) aft 

i V — 

ax Ox 


d (ftv) a (fiv) aift V) f , , 0 0 2 a r ftvi 

« 4. V — _« . ™ v|vt~™ (5 v 2+ 6*^1 +v ’ i 

dt at ax 2 ax 2 ;iv 


dx 

— - V, 
dt 


Loading edge velocity 


For HR >1.0, the advancing loading edge does not influence the method of characteristiu iolution 
in the liquid body. 



^ t 


Figure ~ Advancing leading edge boundary condition. 
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Computer programa were developed and run for both the axial apreadlng 
model preaented above and for a almllar radial liquid apllli apreadlng and 
vaporisation model. Even though aorac coaea have boon run aucceuafully on both 
progrnmRi moat caoea run on theae programa have ahown computational Inotablll'' 
tlea, Thla problem appoara at the free edge boundarlea, and It la apparently 
due to the finite differencing tuchnlquea needed to match the edgo boundary 
conditions to the characteristic aolutlon in the liquid pool. Many different 
approaches have been tried to solve this problem . but none have been successful 
to date. 

Two stably computed runs, with LH^ spilling on a concrete runway, arc 
presented for an axial spread, figure 16, and for a radial spread, figure 17. 
The axial spread resulted from a 900 kg/sec spill of hydrogen occurring from 
an aircraft skidding along the runway at 120 knots. Figure 16 shows the axial 
extent of the spread x measured perpendicular to the spill center line as a 
function of time t at one spot on the runway. The outer curve (square-dots) 
shows the advancing outer edge of the spread out to where It stops and starts 
to recede as the liquid evaporates. The Inner receding edge (circle-dots) 
starts with dry-up near the center at t 9 seconds and then recedes out to 
a final dry-up at x ^ 5 meters and t ~ 11.5 seconds. 

Figure 17 shows the results of a 108 kg radial spill of which occurs 

over a 10 second time span at a fixed position. The center dries up soon 
after the flow ceases, while the final dry-up occurs at r » 4 meters and 
t 17 seconds. For both the radial and axial flow cases, central dry-up may 
or may not occur, depending on the central liquid depth, the liquid inlet flow 
momentum, the rate of vaporization, etc. 

For comparison, the results of simple constant depth slab evaporation 
models are plotted as dashed lines on figures 16 and 17, based upon the work 
of Fay, reference 22. These Fay models assume the leading edge Froude number 
is » 1.4, since V » \^2g6, and hence they overpredict the spread of the lead- 
ing edge, by up to » 20%. Vaporization stops instantaneously with these con- 
stant depth models, and the wider liquid spread yields a faster dry-up time, 
by up to 40%. Based upon these and other comparisons, it appears that the 


63 










Fay models can predict tha extent and time of the liquid spread within the 
limits presented above, With the Fay models, the maxlm^au rate of evaporation 
occurs when the liquid first roaches Its point of raaxlmtun spread, and this 
rate may be up to « 50% greater than the computer model results. 

Because of our inability to make more than a few stable runs with the 
spill, spreading and vaporisation computer programs, it was finally decided 
to use slab models to predict the behavior of the Scenario spills in this 
study. Based upon similar approaches as Fay, reference (22), slab spreading 
models were derived which predict the maximum extent of the spread, the time 
to reach the maximm extent, and the maximum rate of vaporization from the 
liquid spread at this time. These models were developed for an Instantaneous 
axial spill at one location along the spill centerline, and for both an In*- 
stantaneous and a continuous radial spill at a fixed location on the spill 
surface. 

The Instantaneous axial spill corresponds to a spread at a fixed point, 
perpendicular to the spill centerline velocity of the aircraft. The following 
Fay model equations were derived for this Instantaneous axial spill and spread. 


X ■> 


m 

F 

e» M 

1.4 K 

W “ 

m 


• 1 

m_ ■ 
£ 

l-l 

Q. 


•(^) 


*/2 ) 2 /, 


(7) 


where 

■ maximum half axial spread, from spill centerline 

t " time to maximum axial spread 
m 

liig ■ maximum mass rate of vaporization per unit spill length. 
m volume of liquid spill per unit of spill centerline length. 
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The Fay model equations for an Instantaneons radial spill at one location 
on Che spill surface are given below. 



t " 
m 


0.7 


4 r 


m 

3g V 


2 


E 


1/3 


m. 


IT r 


m 




( 8 ) 


where 

r - maximum radius of the spread from Che center, 
m 

0 ■ total volume of the liquid spilled 

^0 

Ag “ maximum rate of vaporization from the liquid spread. 

For a continuous (steady-flow) radial spill at a fixed location, the following 
experessions were obtained for a slab spreading model. 


r 


m 


c 


m 



mjj 


"'m \ “l 


(9) 


where 


r 


m 


a 

Q 


o 


maximum radius of the spread (also the steady-state radius) 
volume rate of the liquid spill 
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The above Fay models for predicting the maximum extent i time, and 
vaporization rate for the various spill typos are most applicable when the 
still rate of liquid fall due to evaporation, Vj^, has a constant value. If 
the term Vg varies with distance and time, as It docs with vaporization due to 
conduction heat transfer from the base surface, then care must be used In 
applying these models. In this case, an average heat flux must be used for 
the extent of the spread at the maximum time, and a trial and error solution 
Is required to match this average heat flux to the assumed liquid vaporization 
rates. Hence, these Fay models only provide an estimate of the liquid spill 
performance, but they do allow comparisons of different liquids as to types 
and scale of the spreading and vaporization. 

For the radial spreads, the instantaneous model was used for spill times 
smaller than the maximum spread time, while the continuous model was used for 
all other cases. In this study, axial spills were only computed with the air- 
craft In motion so only an instantaneous spill model was needed to predict 
this spill performance. 


4.2 Aircraft Fuel Gaseous Dispersion Model 

An aircraft fuel gaseous dispersion model (AFGASDM) was designed to pre- 
dict the position and fuel vapor concentration at any time after the spill is 
on the ground and until the vapor concentration drops below the lowest flam- 
mable concentration. 

Although the literature contains complex grid models using the partial 
differential equatlonr of motion and conservation (reference 23) and traditional 
gausslan dispersion models (reference 24), none of these were suitable to solve 
the problem of aircraft fuel spills. This new, innovative model was developed 
on basic physical principles and knowledge of turbulent entrainment (refer- 
ence 25). The gausslan models are unsuitable because they compute concentra- 
tions based on averaging times of at least 10 minutes, a requirement set by the 
empirical development of the Pasqulll-Gifford diffusion coefficients (refer- 
ence 26). NASA spills of liquid hydrogen Indicated dispersion times of less 
than 1 minute (13). The complex grid models are unsuitable because they re- 
quire implementlon on large computers like the Cray-1 and each run of such 
models costs many thousands of dollars. 
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The approach described here was designed to predict the spatial position 
of the gaseous fuel at each second In time as the liquid fuel vaporizes. The 
concern is the roaxlmvou distance downwind and maximum height above the spill at 
which the gas is flammable. The lowest dctonable concentration is higher than 
the lowest flammable concentration, hence the maximum distance of detonable 
g.as is less than the maximum distance of flammable gns. One second was 
chosen as a convenient interval of time to repeat the computations of energy, 
mass and momentum conservation. The model can easily be adjusted for any 
other interval of time. For the cases reported here, much smaller time 
intervals are expected to increase the computing cost without appreciable 
increase in accuracy of the predictions, while much larger time Intervals 
are expected to reduce the accuracy and resolution of the predictions. The 
accuracy is defined as the closeness with which the model can predict the 
measured results of field experiments. The resolution is the difference in 
predicted results the model can compute for cases with different input 
conditions. 

The model can be run repetitively for each 1 second "puff" of gaseous 
fuel that vaporizes from the liquid. For economy of computer time and cost, 
only the largest puff is followed by the model for the present work. This 
1-second-long puff is largest when the liquid spill has its maximum area and 
the liquid vaporization rate is highest. All other puffs travel smaller dis- 
tances before the gaseous fuel concentration is dispersed below the lowest 
flammable concentration. 

An important aspect of this puff model is that the "observer" travels 
with the puff in a Lagrangian coordinate system whose origin is the center of 
the puff volume. Ambient air is turbulently entrained Into the puff through 
the sides and part of the top and bottom. The behavior of this turbulent 
entrainment is taken from the similar appearing turbulent dynamics of cumulus 
clouds. The gaseous hydrogen dynamics In the NASA White Sands experiments 
were marked with cloud water droplets and observed in movies and in trans- 
parencies supplied by NASA. 


The proportion of puff top and bottom available for entrainment 1 b 
determined by the oacent angle of the puff. There are other puffa before and 
after the largest puff tracked by the model. If the oscent Is vertical, then 
no part of the top or bottom Is an entrainment surface. If the puff motion 
la horizontal along the ground, then part of the top but non«! of the bottom 
la available for entrainment. 

After all the conservation equations are solved, the resulting accelera- 
tion and velocity determine the vertical and horizontal displacement of the 
puff center In a set of Eulerian (fixed) coordinates with the origin at the 
center of the top surface of the liquid spill pool. 

The model is designed to handle four fuels; hydrogen, methane, Jet A 
and JP-4. The Initiation of the model is the maximum liquid fuel area and 
evaporation rate. The evaporation rate Is the depth of liquid lost per unit 
time to vaporization. The Initial temperature of the gaseous hydrogen or 
methane is set at the boiling point at ambient pressure. The model is cur- 
rently designed to accept an input ground elevation of mean sea level. The 
model can easily be adjusted for any other ground elevation. 

The initial puff horizontal geometry is set equal to the liquid spill 
geometry, In Scenarios 2 and 4, the hypothetical aircraft is at rest and the 
liquid spill is circular. In Scenario 3, the aircraft is spilling liquid fuel 
while decelerating from 120 knots, heading into the wind. These Scenario 3 
spills are axial and the model is developed on a unit length of 1 meter down 
the runway. We assume that every other unit length will have similar gaseous 
puffs arising from the liquid. In the axial spills, the puff tracked by the 
model has other puffs in front and in back, as well as on top and underneath. 
This geometry of adjacent puffs is important to modeling the entrainment of 
ambient air because no entrainment is assumed across interfaces with adjacent 
puffs. This assumption is equivalent to assuming the properties of adjacent 
puffs are close enough to those of the tracked puff that mixing across the 
Interface would change no properties being conserved. 

For each Scenario, a model of the liquid spreading and vaporization 
dynamics (see previous Section) is used to provide the liquid evaporation 


69 


race and the apll]. maKltnum radius or maximum axlol wldthi the critical Input 
variables for Che dispersion model. 

Other Input variables are the atmospheric stability • wind speed at 
10 meters • surface temperature and relative humidity. 'Ihe surface temperature 
Is used as the Initial temperature of the Jet A and spills. The relative 
humidity Is assumed constant throughout the layer In which the gaseous fuel 
puffs are entraining ambient air. This layer la no deeper than 762 meters In 
the Scenarios. The wind speed la assumed to have a power law profile of the 
form: 

u(z) - u(z^) (10) 

where 

■ reference height above ground 
u(z^) “ wind speed at reference height 

p ■ power law exponent. 

The reference height is the meteorological standard height of 10 meters. The 
power law exponent Is set equal to 0.2 based on field measurements in a neu- 
trally stable atmosphere. The wind speed at the reference height is set by 
the model user only for a neutral atmosphere. It Is constrained to reasonable 
values selected for the other five stability categories. The six categories 
are labeled A through F. with A being the most unstable and F being the most 
stable (reference 26). The reference height wind speeds in meters per second 
are given in table 17. 

The temperature profile is different for each stability* decreasing with 
height most rapidly for stability category A and increasing most rapidly for 
category F. The stability is chosen by the model user for the mixing layer, 
assumed here to be the lowest 200 meters of the atmosphere* above which the 
stability is assumed to be neutral. 
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TABLE 17. - ATMOSPHERIC STABILITY CONDITIONS 


Stobility Clas! 

Heferoiicfl 
Heiuht 
Wind S|)»i)d 
(m/««c) 

Laiiio 

Hflta (K/kin) 

tettar 

Nuintar 

A 

1 

2 

.40 

B 

2 

3 

-30 

C 

3 

S 

-20 

0 

4 

Any 

-10 

E 

S 

3 

0 

F 

6 

2 

+10 


Tlie puff to be tracked with the model is defined when the liquid spill 
achieves maximum area. For liquid hydrogen or methane, the vaporization 
during the first second determines the vertical thickness of the puff. Tlie 
mass of fuel vaporized is given by; 


"*f " (1-0 sec) 


( 11 ) 


where 

- mass of gaseous fuel in puff (g) 

E - evaporation rate (m/sec) 

Pjg - density of liquid fuel (g/m^) 

A^ - maximum area of liquid (m^) 

The rapid vaporization rate of the cryogenic fuels is our Justification for 
assuming the initial puff is pure gaseous fuel, not yet containing ambient air 
to be turbulently entrained. Hence, the volume of this puff is given by; 

V ■ mf/p„ 

where ® 


V ■ puff volume (m^) 


( 12 ) 



Pg " density of gaseous fuel at 1 atmoapliete pressure 
and boiling point temperature (g/m^) 

and Its thickness Is: 


AZp - V/A^ ( 13 J 

where AZp - vertical thickness of puff (m) . 

For circular spills, the radius is the horizontal measure of the total puff. 

For axial spills, the down runway length of the spill is many times larger 

than the transverse spill width, hence the puff is given a unit length of 
1 meter. 

The initial puff above Jet A or JP-4 is not pure gaseous fuel because the 
saturation vapor pressure of these fuels is small. In order to make the dis- 
persion of gaseous fuel above these conventional fuels comparable with the 
dispersion of the gaseous cryogenic fuels, an Initial puff thickness Is set 
at 7 cm, roughly the thickness of initial 1 second long puffs above the cryo- 
genic liquids. The conventional fuels evaporate into a slab of ambient air 
brought over the liquid at the wind speed for mid puff height for the length 
of time it takes to cross the liquid in the up runway direction. If the re- 
sulting gaseous fuel concentration reaches saturation, then the saturation vapor 
pressure determines the concentration. If the initial concentration is smaller 
than the lowest flammable concentration, then no model run 1s made. 

Once the initial puff is defined in shape, temperature, and gaseous fuel 
concentration, dispersion begins by the entrainment of ambient air and overall 
motion of the puff. The puff of a gaseous cryogenic fuel starts at the boiling 
point of the cryogenic liquid, much colder than its surroundings. The puff of 

conventional gaseous fuel is assumed to always have the same temperature as the 
surroundings. 

Modeling the turbulent entrainment is the critical part. Early computa- 
tions demonstrated that molecular diffusion of heat, mass or momentum was 
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orders of magnitude smaller than turbulent diffusion. This finding was 
supported by the movies of the NASA Wliite Sands hydrogen spills, where the 
vigorous turbulence was marked visually by the motion of condensed water drop" 
lets. The model accounts for two possible contributions to the entrainment 
velocity, the rate at which ambient ulr penetrates the puff-air interface noV" 
mal to the surface. The first contribution Is proportional to the relative 
velocity between the overall puff and the air around it: 


u^ ■ 0.124 


p + (u(=^p)-Up) 


1 0.5 


(14) 


where 

“ first term of entrainment velocity (ra/s) 
w ■ vertical ascent velocity of puff (m/s) 


U(2p) 


U 


horizontal wind speed at height of puff center (m/s) 
horizontal velocity of puff (m/s). 


The constant (0.124) was taken directly from reference 25. In the vertical 
direction the puff is moving through quiescent air, while in the horizontal 
plane the wind is moving air past the puff. The puff is dragged along by 
the wind through the process of entraining air, subject to the conservation 
of horizontal momentum. Horizontal momentum must be conserved because there 
is no net horizontal force applied to the ambient air system surrounding the 
puff. 


The second contribution to the entrainment velocity only applies to the 
cryogenic fuels. It is based on the vigorous mixing Induced just above the 
liquid surface by the high boiling (vaporization) rate. This term is given 
byj 


“2 “ 


(15) 
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• sacoml term of entrainment veloeity (m/s) 

" initial volume of puff (m^) 

E " evaporation rate (m/s) 

This second term was made proportional to tlie evaporation rate. The ratio of 
densities gives the volume expansion of a unit mass of cryogenic liquid as it 
vaporlKes. The puff volume ratio reduces the influence of this term ns the 
puff grows and moves away from the liquid surface. 

The sum of the two terms gives the total entrainment velocity: 


u 

e 


“i “2 


( 16 ) 


where 

- entrainment velocity (m/s). 

The entrainment velocity is multiplied by the time increment (1 second) and 
the entrainment surface of the puff to get the volume of air entrained into 
the puff each second: 


V 

e 


u At A 
e e 


(17) 


where 


At « 

time increment (sec.) 

A - 

e 

entrainment area of puff. 


The entrainment area of the puff is less than the total area of the puff. The 
entrainment area is the sum of the puff side area and the exposed parts of the 
puff top and bottom: 


A 

e 




(18) 
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A ► 
8 


And A, 


J1 


2 

entrelnJug uida araa of puff (m ) 

2 

ftnfjraluing tcp arua of puff (m ) 

2 

ontjralniuis bottom avoa of puff (m ) 


Only if the puff vImon from the ground can tho bottom bo oxpoaod to tutbulont 
untralument of Hurrounding air. The top of tlio puff in partially oxpoaod for 
antrolnineitt oxcopt if it viaua au’uight up. for a c.lroular aitill i 


A " Ai4 (t) 2it r(t) (l‘i) 

u P 

- Ajj - Trr^(t) (Ax(t)/?.r(t)) (20) 

■ Ax(t) 

2 


where 

r(t) ■ puff radius (m) as a function of timu t 

Ax(t) * amount of downwind movement of puff center (m). 

The ratio Ax(t)/2r(t) gives the proportion of the puff top or bottom exposed 
to surrounding air. The other ports of the top or bottom are covered by the 
preceding and following puffs. For an axial spill, there are alinilar relations: 

A - 2Az (t) L (t) (21) 

s p' ' X 




( 22 ) 


where 

L (t) ■ up runway length of axial puff (m) 

X 

L (t) " crosswind (cross runway) width of puff (m) 
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In thp uxlal c|» ,mff Iw, .wlBliborlne puffa m Ua dwuBlnd an<l upwind 

sicicti, na wtill AS Abuve and belpw, 

Snlvlns Ilu. abovu u,„„ap„„ In,.- U,p u,«,„l„,„„, 
oi Air Allpws comiJutaUon of the aiitvained nw«B; 


m 


0 


V p (?! ) 

o'a 


( 23 ) 


whuru 


m 


mass (k) or Air autrulnud duriua a Umo liicrt-ment 


“ danalLy of air at height of puff (g/m^) 

Thin untralnud maan of air la addad to tho oxl.tlng puff „aaa to got tho now 
put nmaa. Tho Initial puff maaa tor tho cryogonlo aplUo waa glvon by 
Muatlon U. Tho initial puff maaa tor tho convantlonal tuela la given by 


m 

po 


VP (z )M 

• _o po P 

R* T (z ) 
a po' 


( 24 ) 


where 

"po " initial puff niaBs (g) 

P - pressure of atmosphere (dyne/cm^) 

*po " initial puff height (m) 

Mp - puff molecular weight (g/niole) 

R ■ universal gas constant (cal/mole K) 

Tg ■ temperature of ambient air (K) 

The putt molacnlar walgbt dopaada on tho molooular troctloaa of gauooua tool 
and air* 


M - M F + M 

P f a 


(1-F) 


(2b) 
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wherp 


II . 


: V 


Mj. • fuel moloculav weight (g/inolp) 

• air molecuJar weight « 28.9 g/mole 
I' " fuel nu)lu I'ractlon 

Tlie fuel mole fraction la an output of the model after the coraputntlona for 
onch time Increment: 


p ^ 

/ ^ ro Mr\ 

I + j._a J 

\ "a/ 


(26) 


where 

^0 “ initial fuel mole fraction 
■ total entrained air maaa, 

Now the information derived from the conservation of mass allows application 
of the conservation of enthalpy. The entrained enthalpy is given by: 


a 

e 


m C T (z ) 
e pa a p 


(27) 


whf re 

“ enthalpy of air entrained during time Increment (calories) 

^pa “ specific heat capacity of air at constant pressure (calories/g K) 
^a^^p^ “ temperature oi air at puff height (K) . 

Conservation of enthalpy proeecds stepwise for each thermodynamic process 
occurring during each time increment. This stepwise solution of the new puff 
temperature between each thermodynamic process provides a higher degree of 
accuracy In specifying the temperature. The theimodynamic processes that are 
handled in the model are: 1) the heating of the gaseous fuel by the warmer 

entrained air> 2) the condensation of water vapor in the entrained air as it 
is cooled by the gaseous fuel, 3) the freezing of the water droplets upon 
further cooling, 4) the melting of these same frozen hydrometeors later in the 
life of the puff, and 5) the evaporation of the water droplets still later. 
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ConcelvflUly the initial cryoganlc tempejratujre (20.3 K) u£ the liquid 
hydrogen can eauHo atmonphorlc oxygon to condonae at 90 K and nitrogen to con- 
donac at 78 K. Xho condoi^aatlou of thcac moleeulea would cause a rapid local 
depression in the pressure and collapse of the volume affected. The heating 
of the cryogenic temperature gaseous fuel to temperatures above these boiling 
points le so rapid, a few seconds, that the latent heat released In the con^ 
dcnsatlon is almost Immediately reabsorbed by the subsequent evaporation. 

There la, hence, no effect on the overall enthalpy conservation. 

Just considering the first thermodynamic process, heating the gaseoun 
fuel and cooling the entrained air, the Intermediate puff temperature is 
given by: 


11 + H 
e P 


T ' ■ 

p C (m 7~nT^ 
PP p e' 


(28) 


where 


T' 

P 


puff temperature resulting from heating of cryogenic gaseous 
fuel (K) 

Hp ■ puff enthalpy before adding air entrained during time increment 
(calories) 

Cpp “ specific heat capacity of puff at constant pressure (calorles/g K) 
The specific heat capacity of the puff Is: 


C -(T )M_r + C M (l~r) 

C (r,T ) - -P^ P f P a a 

pp' p" W 


(29) 


where 

Cpp(l’.Tp) * puff specific heat capacity at constant pressure as a 

function of fuel mole fraction and puff temperature (cal/g K) 

^pfC^p) “ gaseous fuel specific heat capacity at constant pressure as 
a function of puff temperature (cal/g K) 
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In order to compute If water vapor In the entrained air will condenaoi the 
water vapor preasure In the ambient air la found from 


0 - h 10 





}] 


( 30 ) 


where 

e - water vapor pressure In ambient air (millibars) 
h * relative humidity in ambient air 
a « constant " 9.4051 
b ■ constant ■ 2353 

The saturation water vapor pressure In the puff at temperature T’ Is; 

P 

fa - (b/T’)] (31) 

e - 10 L J 

sp 

where 

e__ ■ saturation water vapor pressure in puff (millibars) 
sp 

The entrained water vapor will condense if e < e. If the water vapor 

sp ^ 

condenses, then the latent heat of condensation will provide enthalpy to 
further warm the puff : 


X 

P 


M 


II + II + ra L 
e p w w 


C. (F,l’ ) (m + m ) 
PP P' P e 


(32) 


where 

■ Intermediate puff temperature after condensation (K) 
m^ ■ water vapor mass in entrained air (g) 

* latent heat of water vapor condensation (cal/g) 


79 


Note Che elmllarlty of chle equation with (28), The woCef vapor mass entrained 
and available for condensation Is: 


- e V M /R* T 
w e w a 


where 

" molecular weight of water (g/raole) 

The condensation of entrained water vapor leads to a teat of the Inter- 
mediate puff temperature Tp” in order to see If It is low enough to Induce 
freezing of the water droplets. If 1^" < 253 K, the droplets are assumed to 
freeze. This temperature Is chosen as a compromise between 233 K, when all 
liquid water will freezei and 273 K. when water with many impurities will 
freeze. The presence of ice nuclei in the atmosphere will encourage the 
freezing process, facilitating its occurrence at a higher temperature within 
the range of 233 K to 273 K. There is no need to put this nuclei dependence 
in the model because the required input information does not relate directly 
to the aircraft fuel questions of interest here. Again, the puff temperature 
is recomputed, using the relation: 


T * n 

p 


H + H + lit (L + L ' ) 
e p w w 

C (F,T ) (m + m ) 

PP P P e^ 


(33) 


where 

T^" - intermediate puff temperature after fusion (K) 

^w " latent heat of fusion of water (cal/g), 

The conservation of enthalpy part of the model also keeps track of the melting 
of frozen hydrometeors and the evaporation of condensed water droplets later in 
the dlsperslcm of the puff \.hen it becomes sufficiently warm. The frozen 
hydrometeors ore melted when the temperature of the puff rises above 273 K and 

the condensed water droplets evaporate when e < e . The temperature after 

I ® P 

these processes is T 

p 


Of'trr-"’' ! ■' 

01 1 . 


At the eiid of the conaexvatiot of enthalpy coraputotlons, the final puff 
temperature la used to compute the puff denalty from: 


P 


P 


R#t T "" 

P 


The new volume of the puff Is given by: 


(3A) 


v' 


m (t) + m (t) 


(35) 


The new linear dimensions of: the puff are computed from the new volume assuming 
the previous ratios of the linear dimensions are preserved. For a circular 
spilli 


r 

V' - 
r* ■ 

and A s' ■ 

P 


Az 


Az 

(36) 

P 


n ( r')^ Az' 
P 

(37) 

/ rV' \ 1/3 

(38) 

Az (r'/r) 

(39) 

P 



where 

r' " new puff radius (m) 

Az ' • new puff thickness (m) 

P 

V* “ new puff volume (m ). 
Similarly I for an axial spill • 



(40) 
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or* 

nr *■ r " ‘ , 


p 

h* 

X 

■ ^ 

(41) 

V - 

L' L' Az’ 
X y n 

(42) 

L ' ■ 


(43) 

L' ■ 

y 

L (L’/L ) 

y ' x' X 

(44) 

Az' ■ 
P 

As! (L'/L ) 
p X x' 

(45) 


where the primed variables are the new linear dimensions as a result of the 
puff expansion with entrained air. 

The next part of the model computes the mass motion of the puff. The 
buoyancy acceleration determines both the vertical motion of the puff and Its 
level of turbulence. Telford (reference 25) found that three-fourths of the 
buoyancy potential energy resulted In turbulence In cumulus clouds and the 
remaining one-fourth resulted in vertical motion. Therefore, the net vertical 
acceleration of the puff Is 


A 



(46) 


where 

2 

A * net upward acceleration of puff (m/s ) 

2 

g « acceleration of gravity (m/a ) 

The net buoyant acceleration provides the impulse that must equal the 
change in vertical momentum according to Newton’s Second T.aw of Motion. Hie 
impulse Is the net upward force times the time increment. The initial ver- 
tical momentum If the puff mass times the beginning vertical velocity while 
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die final vertical momentum la the combined maau of puff and entrnlned air 
tlmea the final vertical velocity! 


nipA At 


(m + m ) w' 
P o 


'• m w 
P 


Rearranging, we aolvo for the final vertical velocity: 

m 

w’ « (w + A At) - — 

' m + m 
P « 


( 47 ) 


where 

w' Is final vertical velocity after the net acceleration A has acted 
on the puff during the time increment (m/s). 

The new vertical position of the puff Is found from 

(“ 2 ^) " 

The computation of the new horlKontal position In response to the drag of the 
wind must also take Into account conservation of horizontal momentum because 
there Is no net horizontal force on the puff /air system. Hence, 


where 

u^ “ new horizontal puff velocity after entraining air during time 
increment (m/s) 

and 




((“p ^ “p)'^) 


At 


( 50 ) 
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Thf- ovpi'rtl 1 fiHCPut <u)p1 n of tho puff rul.rttlvu to it-H ntflvtluf? point- 1 h 
uoniputtnl fvoin 


■« tan ^(s<, /x ) (51) 

P P 

Onct! tnc! new vertlcnl position of the puff is computed with (48), the puff 
temperature Is corrected by thi! dry adiabatic. cooliUR it experiences aa the 
ambient pressure drops durluR ascent. The corrected puff temperature is; 

T - T"*' - \'(z - \ (52) 

P P \ P P/ 


where 

Tp « final puff temperature corrected for adiabatic cooling of ascent (K) 

r ■ dry adiabatic lapse rate of the atmosphere (K/m), 

The dty adiabatic lapse rate is used here because any release of latent heat 
from the condensation of water vapor has already been handled in the computa- 
tion of puff temperature during the conservation of eitthalpy. The lapse rate 
of the atmosphere is set in the lowest 200 meter layer for each stability 
class as shown in table 17. 

The lapse rate for class D (neutral stability) Is set equal to tlie dry 
adiabatic value. The equal increments between classes is assumed for con- 
venience, but the sequence is In the correct order with the correct sign, 

4.3 Vapor Purge Model 

This section describes a vapor purge model to be used with Scenario 1» 
small-scale liquid spills into fuselage or wing compartments within the air- 
craft. These small-scale liquid leaks are assumed to be continuous, and they 
may occur from either fuel tanks or fuel lines into adjacent fuselage or wing 
structural compartments. These compartments are assumed to have liquid drain 
holes at their bottoms, so that liquid pooling should not occur In the com- 
partment volume. In order to reduce the hazard of combustion, it is desirable 
to purge the combustible vapors from these compartments at a rate that would 
maintain the vapor concentration below combustible limits. This purge could 
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be aeuompllebod with ram air flow through the various compartinents during 
night, or with pumped air purge If the aircraft Is rm the ground. The 
ob.leetlve (jf this section Is to develop some simple moclelu that will predict 
the required air purge rates au a function of the assumed fuel leakage rates, 

Kor the cryogenic liquids, with small"Scalo leoks Into cotipartinents above 
the. bo.Ulng point temperature, complete vaporlzatlou of the liquid fuel should 
occur at the leakage rate. Except lu the Immediate vicinity of the l..squld 
leak, the purge air rate required to keep the compartment volume below the 
minimum combustible limit of the fuel volume ratio or mole fraction ratio, 
^min' expresseu as a volume flowrate ratio: 



(53) 


where 

V *= volumetric flowrate at atmospheric temperature and pressure. 

Viii “ combustible mole fraction of fuel in the air. 

Subscripts; F refers to fuel and A refers tc air. 

The above equation can be converted to a mass flowrate ratio by multiplying 
by the appropriate molecular weight ratio. 



(54) 


where 

m “ mass flowrate of fuel vapor or air. 

M « molecular weight of fuel vapor or air. 

Values greater than or equal to the above limiting flowrate ratios will ensure 
a noncombuatlble mixture ir. the bulk of the compartment. If the liquid does 
not completely vaporize and some drains from the compartment, then the above 
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atudy raaonraaa, It „aa d.cldad not to attan.pt to modal and aolvo lor air 
purse rculramonta with tho oonvoatlonal tuola. Ganotally ap.„kl„s, porso 
a rtlow radulromonta tor th. lower vapor praaauro oonvontlonal alroratt tuola 
bo loaa than for U.ka ot tho highly volatlla aryog.n,lo tuola. 

4.4 Heat Transfer Analysis in Aircraft Kucl I'iros 

This section is concerned with an aircraft mishap or crash that will 
rcaolt in a lorgo-acalo tuol aplU with combu.tlon coourrlng ad,,aaont to an 
tact tuaolaga oobln. Iho cloaod tusologo la assumod to bo complotoly or 
partially Imboddod In aomo larga-acalo pool hurnlng tlomoa. llte oh,)ootlva of 
s aoctlon la to davolop a hoot tranafer model that could ho uaod to prodlot 
tho tlmo roqulrod olthor to broach tho tuaolago cabio ot to oroato llfo 
t.roatening Intornal cabin condltlono for tho paasongora. 

tatlv'“m?i T',‘ “ «»‘oao„. 

totlvo «d„l Of tho tuaolago cabin of a lo„g-t„„gc aobaonlc aircraft. A» dla- 

ouaaod in soctlc. 3.3.2. tlguro 8 ptoaontod a croaa acetic, . u„d view ol tho 
I LT “r“ T''"‘ “ ‘^Flcal croaa aoctlon of tho 

oomb ai;ir 

bcaod upon thla tlguro S oroaa aooUon. o tUormol modol ot tho tuaolago 
o-w«ll wan dnvolopod to compute tho tranaiont tomporaturoa ot this oroaa 
.action to flame heating ahaotbod on tho outer tuaolago akin. Tho tuaolago 
akin and atrlngora w.ro comblaod „a tho outer t„or,«l node, abaorblng boat 
from tho tlamoa, while tho Internal honaycomb aide panel node okohangod heat 
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by frea eonvectlqn with an air node in the cabin. The flberslaaa Inaulatlon 
waa renroaented by four thornial nodoa in aeriea, with conduction realatance 
between them. Both thermal radiation and free convection heat tranaier waa 
modeled in the air RajJN on each aide of the fiborRlaHH iimulation, An explicit 
aolution of the node temperaturcH aa a function of time waa then Halved on a 
digital computer for a range of abaorbud heat ratow on the fuaelage wkin. 

A cabin window aaaembly and croati auction wan ahown on figure 10. Three 
plastic window panoH are proHont, which should have radiation tranumiaalon 
properties similar to glass. A thermal model of the window panes was not 
developed, however, because it would have required development of the spectral 
wavelength distribution of the various flumes and modeling of the wavelengtVi- 
dependent radiation transmission and absorption of the window panes. The 
task was considered to be beyond study resources and schedule limitations and 
not of the importance to Justify deleting other tasks. 

Detailed thermal modeling of the windows and the flame heat sources would 
yield an estimate of the radiation heat flux transmitted through the windows 
and the time required for the various fuel flames to broach the windows. There 
is an obvious hazard to passengers posed by radiation heat transfer through 
the windows, plus the susceptabillty of the plastic windows to melt long before 
the aluminum structure of the cabin wall. However, a simple solution to the 
problem is to provide heat resistant window shades with a reflective coating 
on the outside surface. In event of a crash in which fire results, passengers 
exposed to high radiation heat flux would thereby have a means of protecting 
themselves. 



5. ANAT<Vf!IS 01'' CUARII fICKNAKTOS 


The vtiHultB of tl)o mialyHCH for four fuel,H nro pryiiuutcd lu-rt* for 
the intenial. lunk (Secnnrlo 1){ the liquid Hplll Hprcudlu):; and vaporization 
which occurti In ScuimrloH 2, ,'} mid A; tlio RUBcoua dlupcrHlou into the atmo- 
sphere associated with Scenarios 2, 3 and 4 and the heut transfer rates and 
the effects of liquid pool fires flurrouudl.ug an Intact fusela^^e. 

5.1 Small Internal Leak Results - Scenario 1 

I'or the amall Internal fuel leaks of Scenario 1, leakajije la assumed to 
occur from the fuel tanks or fuel lines Into fusuluKe or wing compartments 
within the aircraft. The leaks might be due to cracks developing In welds or 
Joints of the fuel system, and they could develop either in flight or on the 
ground. These small fuel leaks will then vaporize in the structural compart- 
ments, and may create a combustible fuel vapor/air mixture throughout the 
compartment. If an ignition source were present, an explosion could result 
that would severely damage the aircraft. 

Running a purge air flow through compartments critical for fuel leaks 
and/or ignition sources is a way to reduce the fire and explosion hazard. 

The purge air flowrate is set so that the fuel vapor/alr mixture Is maintained 
below the lower combustible limit of the fuel. Hence, only a small region 
n.>ar the liquid fuel will be combustible, with an adequate purge air flow 
through the compartment. Normally, fuselage or wing compartments which are 
susceptable to fuel leaks have drains holes at their bottoms to allow any 
nonevuporated liquid fuel to drain from the compartment. 

The liquid fuel spill rates which were established for consideration In 
Scenario 1 for LH^ spills are tabulated below: 


ee 


U - 1.0 k uf^ kR/ri 
l» - 1.0 X 10"^ Rr/h 
11' “ 1 . 0 X 10 ^ kp/« 

Tlu'fU' (miiiR mwn U'RkaRt' wi>vv nlm> uHtui or tlio otlu>r IuoIm (M.'ll , 
im\ .U't A) onuly«o(! iu Mtudy. 

TIir pm’Ro i\U: n„w riH|nlri'mont,n for tlio oliovo fnol IooIwibo vaioK woro 
I'oiimvitKHl u«lt»R tiui purRo rJv nioiio.1 lOVRuufoil in Reotliui .‘i.'l. Thin modivl 
iUiHotncM tliat all t:lio fuel leaUaRo vapovli^iut, ami c.otiiimU'H au a:l)' flow that 
win maintain the fuel vapor mole fraotlon, at: or l^elow the mlnlimmi 

oombuatihle limit In moat of th.e eompnrlmmit volume. The roaultM of th:lf« 
mmlyalH, for the, four fviulM, art' ahown on table IR. wliere thc« air flowrate 
to fuel luakuRo rate ratio in computed to keep the compartment volume in y 
non-combuatlhle atate. 

I'or «mall~Hcale lenkn of tlie boUluR ;!ryoRcnlc fualt, he aaHuiiiption o'^ 
complete fuel evaporation Is Rood, and tb« computed f,ovrntc= v.itlo will h’ 
needed to obtain a iioncombuutlblo mixture, I'or example, for Lli,, fuel 
leakage rate in the range of lir^ to 10'‘2 kg/s would reciuire a tmrge air Uov- 
rate in the range of 3.A5 x to 3,45 kg/a respectively, to obtain a 8,.fe 
gas mixture. For the conventional fuela. the situation va more complex, since 
fuel vaporization cun only occur by diffusion into the air ntream. Th^. licjuld 
fuel may or may not completely vapor It ia a complex prcbluiii to ana- 
lyze, involving liquid areas and flow gf>ometry, vaper and air p^'essates, 
temperature, and vapor mass diffusion :nto the air flow, T!'e table 18 air- 
to-fuel flow ration for the cor./enticMal fueLr uve conservative since it in 
assumed that all of the fuel lerh vaporizes. For either JP-4 -jr Jet A leaking 
at a flowrate of 1.0 x JO’^ Kp/s, an air flcwrnte of about 0.:>y kg/s will 
yield <4 noncc'MbuBtiblc mlx';ure in the compartment. 

For the low vapor in-essure, nonvola/.lle furcc, a certain temperature 
must be present I- , fore a combustible vapor/alr mixtu,e con exist, even with 
no air flow. Pjr Jot A (ketosinc) VAporizlng Into a one. atmosphere air 
pressure, a temperature of >311 J .< (100"F) required for a combusi lble 
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TABJJ' Ifl 


. " AIR n<OW TO FUEI, l.EAKAGK KATJO,S FOR NON-COMHUKTUU.K CAEFF 



vapor/air mixture wich y>0,006. To obtain a combustible mixture with the 
more volatile JP-4 (ganoline) in sea level air, a lower temperature of >237,8 K 
(-32°F) would yield the combustible mole fractions y>0.008. Hence, for the 
less volatile conventional fuels, the temperature and air pressure may deter- 
mine if an air purge flow Is required for a compartment. 

In conclusion, the purge air requirements for a boiling c.ryog c fuel 
leak in an aircraft compartment are more severe than for a less volatile con- 
ventional f «el. The purge air requirements for Ul^ are an order of magnitude 
higher than for the same I,CH^ leakage rate. The conventional fuels evaporate 
by diffusion, and purge air may not be necessary to maintain a noncombustible 
mixture if the fuel vapor pressure Is low enough. Even if the conventional 
fuel leak completely vaporized, the air purge rate is less for the higher 
molecular weight conventional fuels. Hence, the problems of fuel leaks into 
aircraft compartments are more severe with the cryogenic fuels, but they 
should be relatively safe with purge air flows which are reasotmbly obtainable. 
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5.2 Liquid Spill, Spreading, and Vaporisation HeauHa for 
Scenarloa 2, 3 and 4 

5.2.1 Scenario 2 - Radial aptlla . - In the Scenario 2 apllla, liquid flowa 
from either a fuel tank or fuel line onto a fixed poaltion on the bane surface, 
reaultlng in a radial liquid spread. The liquid hydrogen apllla were 

supplied from a single 12,600 kg capacity tank. Three representative Scenario 2 
apill rates were specified, as tabulated below: 

2A - 0.5 kg/s, line leak due to gravity flow. 

2B " 2,3 kg/s, line leak due to pumped flow, 

2C - 900 kg/s, leak due to tank puncture. 

Both the Scenario 2A and 2B leakage rates are small; hence these leaks appear 
continuous with spill times >1.0 hour unless appropriate shut-off valves are 
closed by the flight crew. The Scenario 2C spill occurs over a 14-second time 
span, and this may approach an instantaneous spill depending on the liquid 
spread and vaporization time. 

In this study, the liquid methane (LCH^) tank loadings and spill rates 
were taken as 2.7 times those given for LH 2 above. For the conventional fuels 
(JP-4 and Jet A), the tank loadings and spill rates were fixed at 3.3 times 
those chosen for LH 2 . Hence, the spill times are the same for all of the 
fuels considered in this Scenario, which effectively compares fuel leaks 
based upon the same degree of damage for each of the fuels considered. 

For the Scenario 2 spills, the maximum radius of spread, r ; the earliest 
time of the maximum spread, t^; and the maximum rate of liquid evaporation, 
AgVAP' calculated based upon the Fay models for radial spread pre- 

sented near the end of Section 4,1, Usually, an iterative solution of these 
three parameters is required, to match the liquid evaporation rate to the 
average of a variable heat flux occurring to the liquid pool, In this ease, 
the Scenario 2A and 2B cases were computed using the continuous radial spill 
equations while Scenario 2C was solved using the Instantaneous radial spill 
relations of the Fay model. 
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Two opin RprofuUnr,, mid vnporlzrttlon condltlonn were considered in this 
mid the reumlniiir. Sceiinriofi. The first condlttmi considered was either 
boilliiR vapori station duo to heat transfer from the base surface for the 
cryof»enic liquids, or constant temperature vaporization controlled by vapor 
diffusion into the atmosphere for the conventional fuels. No burnliiR was 
assumed to occur above the liquid piol for this conditions hence, the extent 
of spread and the fuel vaporization rates computed here formed the input data 
for the Vapor Dispersion Models presented in Section 4.2. For this and the 
remaining Scenarios, concrete was chosen as the spill surface, with a sea 
level atmosphere pressure at 297 K (75°F) temperature and a 2 m/s wind con- 
dition. For the cryogenic liquids, a variable heat flux was considered from 
the base surface, while the vaporlzaton rate of the conventional fuels was 
considered constant, based upon a large scale 100 meter liquid spread with a 
2 m/s wind blowing over it. 

The second condition considered was the spill, spreading and vaporization 
of the liquid fuels with pool burning occurring above the liquid spread. With 
ignition at the start of the spill, radiation from the flames would be the 
major source of heat to a large-scale liquid fuel pool. Table 19 shows the 
radiant heat flux qJJ expected to large-scale spills and the rate of liquid 
drop Vj. ’.ue to boiling at this heat flux for the various fuels considered. 

These data are based upon the flame temperatures and the flame emissivlties 
Cj, presented by the A. D. Little, Inc. representatives at the oral presentations 
held at NASA-Lewls on September 2, 1981. For the cryogenic liquids, this ver- 
tically downward radiant heat flux for the equilibrium flame temperature would 
be completely absorbed by the liquid or base surface, and would result in 
vaporization of these boiling liquids. For the conventional fuels, this radi- 
ant heat flux would first increase the liquid temperature, which would also 
transfer heat to the base surface, and finally the boiling temperature would 
be reached. A detailed transient mass and heat flow analysis would be required 
to estimate the time and extent of the conventional fuel spreads with pool 
burning above them. Hence, these calculations were only made for the spreading 
of the cryogenic boiling liquids with pool burning flames above. 
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TAHU-; 19 


. " KXPKCTKD RADIANT IIKAT I-’I.IJX I’ROM VAKIOUH TLAMRS 


Fuel 


‘f 

a ** 

kW/ni 2 

Vt 

in/is . 

LM 2 


0.4 

189.4 

6,0 X 10'^ 

UCH 4 

HiUO. 

0.8 

229.6 

1.1x10'^ 

JP-4 


0.9 

14.7 

2,6 X 10'^ 

Jot A 


0.9 

74.7 

2,6 X 10’^ 


NOTE: 1 Itu tomperBlurc T p and omissivity ‘p of tlmsa flames are based u|ion tho A.D. Little, Inc, con)bustiun 

analysis prGsented at tho NASA-Lewis oral presentations, Sept 2, 1981. 


The results of the spreading and vaporization solutions with no burning 

are shown on table 20. Tor the cryogenic liquids, the radius r and time t 

m m 

to evaporate the hH 2 spills are less than those for LCH^. Vaporization of 
the conventional fuels requires much larger radii and times, although the 
more volatile JP-4 (gasoline) is considerably faster than the less volatile 
Jet A (kerosine) fuel. Scenarios 2A and 2B, which used the continuous radial 
spill model, show the same rate of vaporization as the spill rate m. Sce- 
nario 2C, which uses the instantaneous radial spill equations, shows a lower 
maximum vaporization rate than the spill rate which occurred during the 
14 second spill duration. 

Table 21 shows the results of spreading and vaporization solutions for 
the cryogenic fuels with pool burning flames above. Once again, the extent 
and times of the LH^ spills are less than those "or LCH^, especially for the 
larger rate Scenario 2C spill. Comparing these solutions with flames to 
those without flames, table 20 shows that burning decreases the radius and 
time of the spread, while increasing the maximum rate of vaporization for 
the instantaneous Scenario 2C spill. 

Hence, for the range of radial spills considered here, the cryogenic 
liquids would evaporate much faster, with less radial spread, than the con- 
ventional fuels. Comparing the high rate Scenario 2C spills with no flames. 
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JAIII.K 20. ~ SCl'lNAAtO 2, ItADlAT, SHU., .SP^DINO AND 
VAl'OUI NATION, NO l-'J.AMFlf). 


Scuimria 
No., FiipI 

m 

t (0 

in 

ky/tnc 

Qu 

m 3 

nr/ioc 

m/soc 

2A, LHj 

12,600 

0.6 

117,8 

0.0071 

2.5 X 10 '' 

26, LH 2 

12,600 

2.3 

117.8 

0.0324 

2.5 X 10 ' 

2C, LH 2 

12,600 

900 

117.8 

12.70 

2.5 X 10 ’^ 

2A, LCH 4 

34,398 

1.37 

81.13 

0.0032 

6.0 X 10 '^ 

2B, LCH 4 

34,398 

6.28 

81.13 

0.0148 

4.0 X 10 '^ 

2C, LCH 4 

31,396 

2,457 

81.13 

5.795 

1.1 X 10 '’ 

2A, JP4 

42,210 

1.68 

54.12 

0.0021 

2.0 X 10 '® 

2B, JP4 

42,210 

7.71 

64,12 

0.0099 

2.0 X 10 '’ 

2C.JP4 

42,210 

3,015 

54.12 

3.865 

2 . 0 x 10 ® 

2A, Jet A 

42,210 

1.68 

54.12 

0.0021 

7.0 X 10 ® 

2B, Jet A 

42,210 

7.71 

54.12 

0.0099 

7.0 X 10 ® 

2C, Jet A 

42,210 

3,015 

54.12 

3.865 

7.0 X 10 ® 


i.a O.G 

3.0 2.3 

32 530 

3.6 1.37 

8.0 6.28 

117 560 


237 7.71 

785 100 


the LH 2 would completely vaporize in 32 seconds, compared to an LCH^ vapori- 
zation time of »117 seconds, wliile the maximum hydrogen radius is approxi- 
mately half that for methane. Burning above these cryogenic fuels decreases 
both the radii and times of the spills by comparable ratios for both T.H 2 and 
LCH^. The noburning vaporization times of the conventional fuels may seem 
short, but it should be remembered that these models have a completely flat 
surface, with no small-scale liquid pools. Although we did not solve spread- 
ing of the conventional fuels with flames above, a burning Jet A. Scenario 2C 
spill might last for approximately an hour based on aircraft crash experience 
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TABLK 21, " SCENARIO 2, RADI AT, SPIIJ., SnUiADINfi AND 
VArORI7,ATK)N, RT.AMRfi AlTOVl'T 


Sconarlu 
Nu., Fuel 

m 

kg 

m 

kg/ioc 

Qo 

m3 

ni3/suc 

m/seu 

fm 

ni 

<in 

IOC 

'"evap 

kg/ieu 

2 A, LH2 

12,600 

0 .G 

117.8 

0.0071 

6.0 X 10 ^ 

0.6 

1,0 

0.6 

2 e, LHj 

12,600 

2.3 

117.8 

0.0324 

6,0 X 10 ‘® 

1.3 

1.7 

2.3 

2C, LHj 

12,600 

900 

117.8 

12.70 

6,0 X 10 ’^ 

26 

16 

9U0 

2A. LCH 4 

34,098 

1.37 

81.13 

0.0032 

1.1 x 10‘2 

1.0 

2.6 

1.37 

2B, LCH 4 

34,398 

6.28 

81.13 

0.0148 

1.1 X 10‘3 

2.1 

4.1 

6.28 

2C, LCH 4 

34,398 

2,467 

81.13 

5.796 

1.1 X 10 ‘® 

35 

30 

1,710 




The largest Jet A aircraft spill (the Tenerife crash of two fully tanked 747s 
on March 27, 1977) burned for more than 10 hours. 

5.2.2 Scenario 3 - Axial spills. - For the Scenario 3 spills, liquid flows 
from either a fuel tank or fuel line onto the base surface while the aircraft 
decelerates from a velocity of 61.7 m/s (120 knots) along the spill line. 

The tank loadings and leakage rates are exactly the same as those used for 
the Scenario 2 spills. The axial spread considered was a one-dimensional 
spread perpendicular to tlie spill line. With the high aircraft velocity, 

V^, a constant rate of spill will yield a near Instantaneous spill of 0^ 
(volume per unit length of spill line) at a fixed position along the spill 
line. Neglecting the velocity along the spill line, this instantaneous spill 
will then create a unit width spread, out perpendicular from both sides of 
the spill centerline. In this model, the axial distance, x, will be half tlie 
full spread distance, measured perpendicularly from Che spill center line to 
one edge of the axial spread. Based upon the Fay model equations for 
instantaneous axial spills presented in Section 4.1, the maximum axial Imlf- 
spread distance, x^; the earliest time of the maximum spread, t^; anJ the 
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nmximum rate* of evaporation per unit width of apread, were raU'ulated. 

The anme methoda of aolutlon and lioth the hurnlnp and non-hiirnlnp, eondltloua 
were computed for the Scenurio 3 spllJa, ua waa done for JJeeuario li, Kor 
Seunario 3, only an Inatantancoua axial apill model waa rctiuived, eonaiderint* 
the aircraft velocity. 

The reaulta of the axial spreadlnR and vaporisation solutlona with no 
burnltiR are shown on table 22, while those with pool burnlnR above the liquid 
spread are given in table 23. Tt should be noticed that the low spill rates 
of Scenarios 3A and IB yield very short axial spreads, x^. Considering the 
high aircraft velocity, • 61.7 m/s, the liquid would actually spray upon 
falling and/or striking the spill surface. This would tend to increase the 
extent of the spread and decrease the time to complete evaporation compared 
to those values presented In the tables. This Instantaneous axial spread 
model best simulates the high spill rate cases, such as Scenario 3C, where the 
one-dimensional axial spread extends for a considerable distance. 

Generally speaking, the same observations made for the Scenario 2 spill 
data pertains here. The extent and time of the maximum spreads Increase In 
this order: LCH^, JP-4 and Jet A, the order of decreasing volatility. 

As expected, flames above the liquid spill decrease the time and extent of 
the vaporization, whi-'h occurs at an increased rate. 

In conclusion, a fuel spill from a moving aircraft distributes the liquid 
along a long spill line. The spread perpendicular to this line is then 
shorter in distance and time than for a fixed location radial spill, For 
example, the Scenario 3C, LH^ apill with no flame extends 864 meters down and 
6.2 meters on each side of the spill line, with the aircraft moving 61,7 m/s. 
The total time of spill along the spill line length is 14 seconds, while the 
evaporation time at any point along the spill line is *13 seconds. In com- 
parison, the Scenario 2C, LH^ radial spill with no burning, had a maximum 
radial spread of 35 meters and required «32 seconds to evaporate. Hence, the 
aircraft motion in the axial spills increases the surface area of the spread, 
which substantially decreases the time required for evaporation at any loca- 
tion along the spill line. 
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TARl.M 22. ” AXIAI. SPII.I., SPlUvAniNG AND VAP0H17.ATI.0N , NO V'l.AMKS 


Snuiiario 
Nti., Fuel 

3A, LH^ 

3Q, I.H 2 

3C, LH 2 


12,600 

12,600 

12,600 


3A, LCH 4 34,398 

30, LCH 4 34,398 

3C, LCH 4 34,398 

3A,JN 42,210 

30, JP-4 42,210 

3C, JP-4 42,210 

3A,JetA 42,210 

30, Jet A 42.210 

3C, Jet A 



''e 

Xm 

m/joc 

m 


1.14x10''* 2.6x10'^ 


B.20 X 10''* 
0.2058 


2.6 X 10'' 


2,5 X 10'- 


81.13 0.0939 


2,3x10'^ 10 


54.12 0.0626 


1.68 54.12 3.48 x 10® 7.0x10'® 2.3 217 

7.71 54.12 1.60 x 10'^ 7.0x10'® 6.6 400 


7.0x10'“ 203 4,330 


Initial aircraft velocity along the sgill line, ^ 61.73 m/sec (120 knots) 


"’tVAP 

ky/msuc 


1.37 81.13 6,22 x 10'® 1.0x10'® 0.06 0.8 0.064 

6.28 81.13 2.40 x 10 ''* 1 , 0 x 10 '® 0.16 I 1.6 0,134 


41 1.85 


1.68 54.12 3.48 x 10'® 2.0x10'® 0.6 29 0.0018 

7.71 54,12 1.60 x 10'^ 2.0x10'® 1.5 53 0,0046 


2.0x10'“ 53 680 0.165 


0.00025 

0.00062 

0.0221 


5.2.3 Scenario 4 "• Radial spills, catastrophic crash . - The Scenario 4 
spillu are the result of a nonsurvlvable crash of the aircraft, In which all 
of the tanks rupture and spill instantaneously, resulting In an assumed radial 
spreading and vaporization of the liquid fuel. Scenario 4A represents a 
crash shortly after takeoff, where the nearly full tank loadings would be 
spilled instantaneously, while Scenario AB simulates a crash Just before 
landing, with an instantaneous spill of *10% of the total tank capacity. 
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TAHi.K 21 . - AKIAI- RNU., Sl'KIiAni Nf! AND VAl'OKI 2 ATI i>N , I’l. AMI'S AllOV)': 


Scoiiiirio 
No., 1 mil 

Ml 

k|| 

m 

kg/stiu 

in» 

a<i 

in3/iii 

Vf 

ni/si'i: 

" III 
III 

*111 

Sin: 

"'i;vAi' 

ku/mmii: 

3A, LHj 

12 ,GUQ 

0 .G 

117,8 

1.14 X lU '* 

G.D X 

10-* 

O.QO 

0.4 

(I.U42 

3Q, Llfg 

12.GQ0 

2.3 

11 '/.8 

0 . 20 x 10 '* 

G.Ox 

10 ^ 

0.12 

11,7 

0.104 

3C, LHj 

12,600 

900 

117.6 

0.2UOG 

0.0 X 

m 3 

4,4 

7.7 

3.74 

3A, LCH 4 

34,396 

1.37 

61.13 

0.22 X 10 ** 

1.1 X 

10'3 

O.UG 

0.6 

0.000 

3B, LCH 4 

34,396 

G.28 

81.13 

2.40 X lO '* 

1.1 X 

103 

0.10 

1,4 

0.137 

3C, LCH 4 

34,396 

2,407 

81.13 

0.0939 

1.1 X 

IQ -3 

0.0 

10. 

4.94 

Initial aircraft valocity alony tliu siiill line, 

~ G1.73 in/soc (120 knuts). 






The total fuel Hplllago. for 4A for LIl^ was 21,600 kg. and the aatne weight 
ratios as used in Scenarios 2 and 3 were used here for the various fuels. 

Since all these spills are instantaneous, the Pay model eciuations for 

Instantaneous radial spills presented near the end of Section A.l were used 

in this analysis. The maximum radius of tlie spread, r j ‘he earliest time 
of the maximum spread, t^; and the maximum rate of vaporization, were 

computed for both the conditions of no flames and pool burning flames above 
the liquid spreads. These data are presented in table 2A for the no-burning 
condition and in table 25 for the liquid pool burnliij c ndltion. 

Once again, these data show the same ranking of the fuels relative to 
the radial extent of the spread and time for vaporization as occurred for 
the Scenario 2 radial spills. „ith no flames, the maximum Scenario 4A spills 

would find a maximum radius of 49 meters and a vaporization time of 19 secondn 

for UI 2 , compared to 77 meters and 140 seconds for l.CII^. With flames above 
the liquid spill, a 37-meter maximum radlous and a 22 second vaporization time 
is predicted for LII^, compared to 43 meters and 43 seci>nds for hnt^, when the 
maximum Scenario 4A spill is considered. The extent of spread and time of 
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TAIH.K V\, SCKNAUXO A, ItADIAI. SPII.L, MJ’WKADiNC AND 
VATOUIJ'ATXON, no FJ,AMKfi 


ScQiiarlQ 
No., FkqI 

m 

Ku 

in 

kfi/nee 

Un 

m3 

Uo 

ni3/ifli: 

m/ioc 

U\\ 

III 

■ — 

’ill 

101 : 

• ■— 

"'liVAl' 

ko.iiii; 

4A, LHg 

21,000 


304.74 

• 4't 

1,9 K 10'^ 

49 

39 

1,040 

4Q, LHj 

2,100 

4IJ 

30.474 

i|> 

2.6 X 10’* 

20 

19 

210 

4A. LCH 4 

00,908 

ifti 

139.08 

.U 

1.0 X lO '* 

77 

140 

790 

4B, LCH 4 

0,890 

(jg 

13.808 

tA> 

1.0 X lO’*^ 

29 

04 

170 

4A, JP-4 

/2.300 

(XJ 

92.77 


2.0 x 10® 

170 

890 

100 

4B, JP-4 

7,230 

lAJ 

9.277 


2.0 X 10® 

74 

600 

20.8 

4A, Jot A 

72,360 

00 

92.77 

00 

7.0 X 10® 

406 

4,790 

28.2 

4B, Jet A 

7,236 

CD 

9.277 

00 

7.0 X 10® 

171 



2,690 

0.02 


TABLE 25. - SCENARIO 4, RADIAL SPILL, SPREADING AND 

VAl»ORIZATION, FLAMES ABOVE 


Scenario 
No., Fuel 

m 

kg 

ni 

kg/soc 


60 
in 3/sec 

m/sec 

rni 

m 

soc 

i^EVAP 

kg/iec 

4A, LH 2 

21,600 

9i 

304.74 

ijg 

6,0 X 10 ® 

37 

22 

1,830 

4A, LHj 

2,160 

gg 

30.474 

ou 

6.0 X 10® 

16 

12 

330 

4A, LCH 4 

08,968 

gu 

139.08 

ou 

1.1 X 10® 

43 

43 

2,660 

48, LCH 4 

0,896 

OJ 

13.908 

ou 

1.1 X 10 ® 

18 

24 

466 
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vnp<irl/.af lull fur i lu' Sruiuirlo /\h cunvc'ut Iniml fn(>] HplUh ^.j-yator than 
iliudu abnvi' for tho no flniiio cnmlitlnn, tahlo ;>/i , and for (ho flame condition, 
wlioro i^xpor louoo ahowa hum timofi proator tliaii an hour an ropronoutat Ivo of 
Si'oimrlo hh op 11 1 m of Jut A, 

Tho ahovo prod let Uuih of .Seonarlon 7. tlimnp.h a apUlH have ahowii that 
tho oxtont of Huroad and tliiio of vapor I Knt Imi la inueh Iomm for the eryof'ou'le 
luol.fi than the eouvont lonal fuolfi, In hotli tlio hurnlnp, find no hurninp, coii" 
ditlon. Comparing, tlio eryoRcnie fiu'lti, the oxtont of the upro/id for 1,11^, can 
hi) down to half an far, wlill.o tho ovaporntlou tlmu can ho ono"fourth to ono" 
half thfit of U;il^, llone.o, from the atandpoint of minimum liquid spruad and 
tafjt ovaporatlon, liquid liydroRon lu the nafest fuel, t>reKcntinR leant hazard 
to t!ie aurroundini'a In event of a craali. 

5.3 CaseouB hiapursion Rusulta tor Scenurloa 2, 3 and A 

A complete llatlnp, of the paBcouH dlsperalon renultH (for LIL, LCU, . and 

fueln) for ScenarloH 2, 3, and A Is tjiven in Appendix A. A auimnury cf the 
retmlts 1« presented In table 26. The Hlgniflcanee of the results are dls~ 
cusaed in this Section. 

Scenario 2A is a small circular spill. The final time, downwind distance, 

and helfiht above ground are listed for the different fuels In table 26 for 

comparison. The listing for the LH 2 case in Appendix A shows that both the 

buoyancy acceleration (column 5), and vertical speed of tho puff rise from 

zero to maximum and then decrease. The overall anple of ascent is about 
o 

.30 . The first entrainment .speed (Uj^ In column 10) also increases from 
zero to a maximum before decreasing, following its dependence on the relative 
velocity between the puff and the surrounding air. Thu second entrainment 
speed (u^ in column 11) decreases monoLonicnlly after its initial maximum 
because the volume dilution quickly reduces the creation of turbulence bv the 
vigorous evaporation rate. The puff radius and thickness increase mono- 
tonlcally as the puff continually entrains air. 

The next two sheets In Appendix A, a plot mid listing for the methane 
case for Scenario 2A, show the LCll^ puff lasts almost twice as long as the 
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hytlvoKf'ii faiK.i and iravnln alniiR the K»’<'iuid f(»r iiumt nl' Uh Ult't liiu.*, An 
oiitVii liiiiioiU' coul^^unnti, ihn Inifirtl nnpni Ivt' Itnoyancy lifronuat ii] Ij'liiJv inu’iiivi* 
and UflH Uif puff n 1 t'w inntt^rn r>ff tlin f^rmmd. 'I’hn 1 Inf lap. nli'»wn now 
onlvailniutMil; dona not ot't'ur tliroupli tin? liottoiii of tin? puff wlton It travolo 
ala)np, the pronnd. 'J'lit* pnnoi’nl, heliav;l<n' of aai'ondln)' liydvopon pnfla /md proniiii- 
it>vo.l, travt’I of motliano puffa i>rod;lott>d liy tho tnodol :1a anpportod by tlio 
liquid bydropcin tUHtH conduutod at Wliitu Sandu (13) and tba liquid natural 
Raa tuHlH conducted at China Lake (14). 

The next plot and llatlng Appendix A allow how the .lP-4 In ('CCii<> '*o 2/v 
never rlaea off the ground because of its negative buoyancy. It travel'} 

Hlowly for a long time because It huga the ground where the wind apecd in 
very low. Tn fact, the power law wind profile requires that tb •. rind rpu<?d 
go to zero right at the Hurface, a boundary condition .■r.poaet? In all 

fluid flow models that have friction at the walls, r.ero is also no contri-" 
bution to entrainment from the second term because it Is computed only it 
there is a vigorous evaporation rate. There is no mixing of air throngh tha 
bottom surface becaus^i the JP-4 hugs the ground in stable air. The second 
line for time equal to 1 second in the Appendix A listing for J)’~4 shows a 
slight decrease in the puff radius and half thickness before these puff 
dimensions monotonlcally increase. The initial decrease is not real, but 
results from computational error in the program. Tills error is small enough 

V 

to neglect in assessing the model results. 

There are no plots or listings for any Jet A cases because no flammable 
concentration is ever produced when the surface temperature is 294 K (70°1'). 
Only if the surface temperature is increased to at least 314 to 322 K (105 - 
120®?) can a flanmwble concentration be produced from the increased evapora- 
tion rate. 

The results for Scenario 2B (for hydrogen) in Appendix A show how a 
bigger spill inoreases the trajectory time, dimensions (also see table 2fi), 
and ascent angle. The JP-A puff for Scenario 2H lasts 16 minutes as it bug;, 
the ground and slowly moves almost 400 meters. 
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Thfi Scenario 2C for hydrogen ahowe the puff moving 575 inetoro above 
ground during a devniwlnd travel of 411 metora. In cniiiparlaon, the methane 
In Scenario 2C travelo along the ground 713 metern and doonn't get off the 
ground except briefly ,iuet before It dilutee to lt« loweot flammable concen- 
ratloui In the Hating for this case In Appendix A the aucent angle of 05“ 
at the one second time results from the vertical expansion of the puff before 
It starts moving downwind* 

The Jl*"4 case for Scenario 2C shows that vapor from u large (143 m radius) 
spill of conventional fuel can theoretically move u long distance (6.8 km) 
downwind close to the ground fer u long time (over 3 hours). In reality, the 
JV-4 vapor would disperse more quickly when it moved off the downwind end of 
the runway, which probably won't be longer than about 4 km (13,126 feet). 

On an airport the trajectory of the flammable gaseous fuels in all cases 
would move towards Incoming aircraft, Safety considerations would dictate 
that measures must be taken to prevent any aircraft from flying, or trying to 
land, through a cloud (many puffs) of any flammable gaseous fuel. 

For the axial geometry Scenarios 3A, 3B, and 3C, the model follows a puff 
emanating from a spill of 1 meter unit axial length. The spill rates in 
Scenarios 3A and 3B are small enough to allow an aircraft to spill fuel over 
the entire length of runway. This might amount to 4000 unit spills, each 
starting with a 1-meter length, following each other downwind. In addition, 
for each unit length of spill there will be puffs that precede and follow 
this largest puff being tracked by the model. 

Data for the hydrogen case for Scenario 4A in Appendix A shows the 
flammable puff moving off the ground in 20 seconds after spreading to 54,2 ni 
radius. It persists for 173 sec. and during that time rises to 782 meters 
above ground before being diluted tc less than 4 percent hydrogen. The rapid 
rise rate means the vapor cloud from a hydrogen spill will not endanger 
persona or property In the surrounding area in event of a catastrophic crash 
of a LH 2 fueled airplane. On the other hand, the methane in Scenario 4A 
moves along the ground for more than one-half hour to cover 851 meters, and 
it only rises off the ground daring the last 1 to 2 minutes. The Jl’-4 gas 
puff could theoretically move 10.4 km downwind along the ground In almost 
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5 l,„ur„. lUwvar. aH.r „ to, «,„t. U would „,„vo on t„„ oud of t|,„ vuuway 
o».d 8r„uu. wouu'. uto. Thuvefor., „,o bo„„Uary uunUlUouo will „,,„,„bry 

uhuuBO l„vu11d«l„B 11, on., ,U„t„.u.o„ ,„„1 t,.„vol ,l„o. 

llowovur, tl,o potoutlnl o»l„,a u.o 

()):opi>K'i;y onto thcj .luiul ot otlicv uwuux’h, 

noto of tI,o movcnout and nlzo of tl.o putfo for oooh of tl,„ f„ol« „ro 
l.,-o„o„u.d 1„ flguroa 18 and 19. pinuru i« „,,owo nUltudo vorouu drift dln- 
tanoo with Indications of olapsud tlm, ,i,nrkod on tho path of each of the 

fuels. Fleurc 19 lllustratos the prowth of the puff diameter as n function 
of time for all three fuels. 

The hydrogen case for Scenario 4A was also run with each atmospheric 
stability class. As shown in table 27, from data in Appendix A, the puff 
rose least and traveled furthest downwind in the least time with stability 
classes 3 and A (C and D) . The wind speed was at the maximum 5 m/s for both 
classes and 5ts effect probably outweighed the effect of temperature profile, 
determined by the lapse rate. 

Ihe effect of relative humidity on the Scenario AA LM„ spill can be 
«oen in the bottom part of table 27. More water vapor reduces the travel 
time and downwind distance monotonically but the maximum height is achieved 
by a puff rising through air with 50 % relative humidity. 

5. A Heat Transfer to the Passenger Cabin 

Consider an aircraft crash with a fuel spill beneath and flames around 
the fuselage, but witli an intact and closed passenger cabin. For the various 
iLiels considered in this study, what time period would he required for the 
flames to breach the cabin, or create a life threatening environment for the 
pnssongers in the cabin? In Section A.A. a thermal model was described to 
compute tho transient temperatures through the fuselage cabin wall cross 
Hcetion. The oblectlve of tills section is to present results of an analysis 
ot tills cabin wall model to estimate the survivability time for passengers 
inside the fuselage cabin imbedded In the flames of various fuels. 

In order to compare the various fuels, the heat absorbed on the fuselage 
sU^ln from the adjacent flames must be estimated. A large scale fuel spin 





Altitude (m) 





t (sec) 


Figure 19. Flammable puff diameter va time for liquid hydrc ^en, liquid 
methane, and JP-4 (Scenario 4A, wind speed e 2 m/s). 





TABLE 27. - GASEOUS WSPERSION RESULTS FOR DIFFERENT 

ATMOSPHERIC STABILITIES AND RELATIVE HUMiniTIES 


Fufll " Hydrpgen. Scenario 4A 


Stebllitv 

6 

6 

4 

3 

2 

1 

6 

6 

6 

6 

6 


Relative 

Humidity 

(%) 


Finsl 

Time 

(tecondt) 

Powiiwind 

Piitance 

(meters) 

Height 
Above Pro 
(meters) 

173 

612 

782 

148 

680 

537 

119 

703 

347 

119 

702 

344 

153 

635 

557 

186 

536 

739 

267 

842 

763 

222 

684 

774 

173 

512 

782 

161 

431 

V70 

138 

357 

718 


was assumed which would result In a large scale pool burning. The fuselage 
was assumed to be directly Imbedded in tl . flames, so that both radiation 
and convection heat transfer occur to the outer surface of the fuselage 
skin. The incident radiation heat flux from each of the large scale fuel 
flames can be taken from table 19. The radiation view factor from the 
vertical fuselage skin wall to the large scale flames should be -0.5. Al- 
though the bare aluminum skin could have a low emisalvlty, - 0. 1 , a painted 
skin or one subjected to corrosive flames could have emloslvlty e- 0.9, which 
was assumed in this analysis. Hence, the absorbed radiant heat flux on the 
fuselage skin q*^, is shown on table 28 for the various fuels. 

With flames next to the external fuselage skin, turbulent convection 
heat transfer will occur from the hot flames. The turbulent convection heat 



TAHI.K 28. AllSOUmil) I'USELAnii SKIN HEAT FT.UX IN THE VAUIOUS IT, AMES 


fuel 

qj; 

“R 

kW/m^ 

a" 

kW/ni^ 

q" 

kW/m^ 

LH2 

180.4 

85.2 

31.8 

117.0 

LCH^ 

229.G 

103.3 

27.2 

130.5 

JP-4 

74.7 

33.6 

18.2 

51.8 

JET A 

74.7 

33.6 

18.2 

51.8 


transfer coefficient, h, opposite the vertical fuselage skin wall, is estimated 
as hss23 W/m K (4,0 Htu/hr ft^ °R) . Using this heat transfer coefficient, 
the convective heat flux, q^, to the skin is shown on table 28 for each of 
the fuels. The estimated maximum heat flux absorbed by the skin, q", was 
then obtained by adding the absorbed radiant and convective heat fluxes, which 
yields the data in the last column of table 28. 

The transient thermal model of the fuselage cabin wall, briefly described 
in Section 4.4, was run over a range of absorbed heat fluxes, q", on the 
fuselage skin. The solid curves of figure 20 show the outer fuselage wall 
temperature as a function of the absorbed heat fluxes labeled on those curves. 
Plotted as dashed lines on the same graph is the temperature response to the 
estimated maximum absorbed heat flux q" for each of the fuel flames considered 
in table 28. The dotted curve at the bottom of this graph shows the Inner 
cabin wall temperature, which only rises one or two °C (two to four °E) over 
this short time span. 

When the aluminum skin and stringer temperatures are in the range of 
430 C (805 F) to 480*^C (896°F), the fuselage should collapse due to its own 
weight. This fuselage collapse temperature band is shown in figure 20. Based 
upon the maximum expected absorbed heat fluxes from the various flames, the 
fuselage should collapse after v. AO seconds in a methane flame, and after 
~50 seconds in hydrogen flame, and after -120 seconds In a JP-4 or Jet A 
flame. The Inner cabin wall temperature barely rises over short time spans, 
and hence this temperature is no threat to human survival. 
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To ovftlufito tlif* poHHlIiTHty of furtolrtRO eoXlapHo In larRO noalo h;i 111 
f lniiiPfi, t;lio oxprototl Imrnlnp, timco for flo.cnnrio ?,(1, !\ti and AH rndlal up 1 11a 
wldi flniiK'ii nlmvp, /ui p.lvon In tnbloa 21 and 2f) wero lnvi-*Ht1p;ai:Pd. For M!^, 

Ulio fuHC'lnpo pollapao tliiio 1", " .'iO socondfi, timch Rroator than tlio larp.oHl 
profUc.ted burn timo of 22 Hoconda. For tha fuHolnnc oollapHO timo la 

->-40 Hoooudfj, which is comparable to the l.arp.eHf predicted burn time of 
43 seconda, For either JI’-4 or .TliT A, tlie fuselape col.IapHe time of 
~ 120 seconds la much less than the expected conventional, fuels burn times 
of approximately an hour. Hence, based upon this study of larpe scale fires, 
the fuselage structure should survive a hydrogen fire, be marginal for col“ 
lapse In a methane fire, and should definitely collapse in J!’-4 or JFT A 
fire about the fuselage. 

Although a thermal analysis was not done on the cabin windows, a few 
observations should be made. A good fraction of the high radiant heat flux 

2 

from the flames would be transferred through the plastic windows, say ~ 50 kW/m . 

2 

Since a radiant heat flux of ~10 kW/m can cause skin burn within a few seconds 
exposure, the passengers would appear to be at risk to radiation burns through 
the windows, especially with the cryogenic fuel fires. Although the thernuil 
capacity of window panes appears to be similar to that of the fuselage wall 
cross-section and only a fraction of the radiant heating will be absorbed by 
the window panes, the softening temperature of these plastic window panes 
should be in the range of 90“C (193“F) to 150“C (302°F). This low structural 
failure temperature of the window pane materials may mean that flames would 
breach the window panes before the fuselage would collapse. Hence, because of 
the possibilities of radiation skin burns through the window panes and struc- 
tural collapse due to the low softening temperature of the plastic window 
panes I a transient thermal and radiation analysis should be done on the cabin 
windows and means of alleviating the high heat flux (such as reflective coat- 
ings on the outside surfaces of the window panes) Investigated. 
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(JONCLUSiONS 


Com;luoiouu ruacUod aa a vuwu.1l o! tliiw uvawli Uvu haaavd atudy ol vavlouw 
aivcval't I’uula uve aw lolloww: 

1 . Cryoguiiic £uelw will bu cavviud lii I'uwulagu tuukw which iiihuvuntiy 
are more resistant to dumugu than is thu rest of the fuselage. 

2 . I'uselage tanks for cryogenic fuels are less susceptible to damage 
which would result in spilled fuel than are tanks for conventional 
fuels which are located in the wing. 

3. Tanks of conventionally-fueled aircraft are vented to the atmosphere; 
therefore a combustible mixture can exist inside the tanks. Cryogenic 
fuel tanks are pressurized with gaseous fuel. No oxygen is present 
within the tanks so there is no fire or explosion hazard. This 
represents added safety for cryogenic fuels. 

The extent of spreading of the liquid fuels, and the time required for 
them to vaporize, varies Inversely with the volatility of the various 
fuels. Ranking the fuels In the order of increasing vaporization time 
and extent of spread, for conditions both with and without flames above 
the pools, yields the following orders LH^, LCH^, JP-4, JET A. 

5. The probability of fire in event of a survivable crash whure fuel is 
spilled is quite high (greater than 80%), regardless of the type of 
fuel involved. 

6 . Ihc ground level hazard to the surrounding community represented by 
spilled hydrogen in event of a crash where fire does not inmiediately 
result is dissipated within about 40 seconds in the worst case 
analyzed. Since a flammable cloud rises to a height of about 800 m, 
it poses a threat to low flying aircraft for a few minutes. With 
LCH 4 the ground level hazard would persist for more than half an hour 
and the flammable cloud would be blown 851 m before it is dissipated. 

It therefore poses a substartial threat to people, vehicles, and 
structures on the ground ir this area. Tlie cloud of flammable vapor 
from the corresponding JP-4 spill never leaves the ground, persists 
for 4.8 hours, and theoretically could drift 10.4 km. Jet A and JP-4 
vaporize by diffusion at the normal ambient temperature condition 
analyzed; however, Jet A does not have a high enough vapor pressure 
to form a combustible mixture at standard atmospheric conditions. 
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/. A flrt* I’tsuiU (-tiiK from (’raoli of n Ul./ fm'loil olrora) t. wi ll lie vary 
mall III aroa ami aliort In dnr ' on'Voiiipariol t:o otlior fnolii. 
I'aaaoiif^o.ni In Ml,^ alrorafi can luohalily luoit mirvlv** a cranlt fire hy 
Ht-ayliiK In tliof r' acatfi until the fire hnnin out (lean than 10 nccomlfi 
for the womt eaae Htmilmi). With .all tlie other fneln a fire woiilil 
vofinlre Immediate evae.nation of the ftiHel.ane, 

H, hlqnld drain holea and inirjje air flow will rediie.e the hazard of flren 
eauaed hy wiia.ll fuel leakn Into aircraft coiiipartmenLii for hotli tliu 
eryo{jeiile. and conventional fueln. 

y. Conveutlonul gausaian diapersion modeln cannot be uaed to predict the 
time/spacu dynamics of gaseous hydrogen and methane orlglnatiiij; from 
cryogenic liquids which disperse sometimes in l__a than one minute. 
Complex, expensive models using tlie partial differential equations of 
motion for a large number of 3-dlmcnsionul grid cells are unncces" 
surlly detailed and expensive, Thu methodology of a lumie-lcal model 
developed for this *audy to solve the energy, mass and momentum con- 
servation equations for individual puffs of gaseous fuel evaporated 
from the liquid spill has proved to be very satisfactory, liacli puff 
is created during a time span of one second and it Is followed at 
every one second Interval until it Is no longer flammable. l-‘or the 
subject study the model was run on the puff produced by the maximum 
evaporation rate from the spll.l wlieti it has its muxlmum area. Any 
other puff can be followed. The model accounts for all thermodynamic 
processes of water vapor phase changes and sensible heat balance. 
Molecular diffusion is negligible compared to turbulent diffusion. 
Turbulent entrainment Is modeled on the current understanding of the 
similar turbulent entrainment occuring in cumulus clouds. 

10. Considering the spill, spreading, vaporization, dispersion, and burn- 
Ing phenomena of these fuels in the subject crash scenarios, liquid 
hydrogen (LIL) fuel appears to offer many advantages when compared to 
LCll^, Jl>~A, or J12T A. 
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7. Ki:COMiV|,N))ATJONti KOU rUTUKK WOKK 


TUii ).'tH;umiiiuiul,iLi(iH!i j'or 
pusuibi.i! uXrci'aJu lui-lu uan b 
and I ueuLii. 


liitiuc! Work Jit ti vt! Uo Itiiiiiivd 
0 dJviduil iuLo LWi) catu'yto'iioj ; 


iU)iiipaj’.luoaii ol 
aitalyUcaJ auuUoa 


1 . 


2 . 


3. 


4. 


5. 


6 . 


/.J Aiuii.ytit’uJ, SLuUiua 

tliL aixn, apioadxiiy and vaporiauUJtiii pry};ruins, TIh'mc Drot-rim- win 
ovlJo moic ocoural-., foi.voaonlavfvo data lot tho I'ala,, aud 
lotloimaiioc l)rcdlctiou al tlia ooaaaijuoiiooa ol liquid apiUt:, 

lucludo lioat utaiialat duo to aurlaoo conduction, llauio radlitjon ml 

c“iu": ::;-tS:'"n!,’!;‘if uiuita edacity 

(-Lict-Ub 111 Llio litpitd vapuriiiation modol. '' 

aru‘'LJ^aLni?L“an'tiv “P-^^dang model 

may sltiixliuaialy impedt. u cryogenic liquid tiprcadiiig IJow. 

H^auous dispersion model with input conditions nrovidod 

UJJ, Arthur D. Little Co. spills ol: liquid hydrogen in 1958 
( eteience 12), and the LNG spill test program conducted by Lawrence 

= dSunf 

Run the gaseous dispersion model over dJifereiit Hurtne.. ..... 
vatylns ftom 233K(-40-l) to 322K (120-« and nutlacu iLlalL'l va^^ 
Ills Iron, sua level up to 1.5 km («922 ft), and wind apeoda vaivl.m 
from dead calm to 30 m/s (approximately 60 inph) . ^ ^ 

Run the gaseous dispersion model for spills of fuel irom storane 
ground tojogrlphj? ^P^^citic surrounding 
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h..| Liiid t )(i' dl(i|u rMluu imulc I (-u -i) lu j iluu;; nj 

r.iitjauiiu'Ufi lilu! dil.,riiu‘, hy.lj-nj>'ii ililml.lr and n|,rml.al waii.n. 


I 

f'.ii 


a. 


I't'i-i arm aiM lyUiia I atndU-;. al ihr .'/iin.tad 
Mania aiiiitihi vJt leu aJ dUli-ti'iil I iia I llama;, 
avaluiiaa l l.itiu; damar.a La ilia .iJi. i.iil i.Lnii 
ill ill! I'xita,- Lad t raali ta'aiiar la, 


I Kuna I rlii|)al a I ill a; 

al v.ii iuiia ;,aala 
'Ilia al ran V i S' i nr, 


. and 
I i a : . I I 
liilraa'lir.i 


y. 


DavuJa|) ail ana I y t. lea I iiiadal La |ii'adJi'l 
Mu:; and Llin u,Muniamt. laiiiiiai ai nran al 
uvalsmLu tin.' liaKai'iln La pantiani’am Irani 


Ilia I ran iiiiii (ad ladlani liaal 
I 111' i-alrlii wliidaw i»aiian i a 
I laiiii'n aiilnldi' till' liii.alar,a. 


1 ,. 


7.2 l';xi>ai'i.iin'iil;i.l TaaLn 

Moi’u 111704 * 1111 X 1 ; inuJ. apllj ts'iiLii, boLli Iniriruij’ and iiun-)HiriiJiir,i waiild 
pi.ovido valuablu uximrioiiuo and dat.a i.a siampaia' arsalui L analyLli'al 
moauJ,.y ot dm wplll HpjHiadlng, vapurJxntlmi mid dinporsiim plmiionmna 
J.bsi analyLJ.cal iiiuUo.l.!i can be uianl La lml.p ski!ii}.;n Lbs; LsniLii. Jn 
udUiLion, but’iiing Lcata could proviilsi uxpui.inmiaaJ vuriMcaUan al 
Uanic Lciiipeniuit’os and radiaLion Muxeii 1 ‘rotn varionn hIus; couM ar,i..'r ■ 
Lionti, ibu opill Loatti should Inc Jude boLli large "iunLantancuun" 
spills, and upillu to siinulutL pIpcUiui lupturo sdmiv lung duraMunn 
arc involved. 


2. Some ilow triction tests oi; a boiling cryogenic liquid spreadinr, on 
a hot suriace should be run to see il the boiling ))hcnomena signii- 
icantly impedes the liquid spread, 

3 . Testa oi (dm radiant heat transmission, tetuperature rise, and 
stiucturul iailurc temperatures should be run on cabin window panes 
and typical aircraft fuselage structure to evaimae the passenger 
hazards, especially with cryogenic fuel flames adjacent to the 
lubelage. 

Crash tests of surplus aircraft with represenla 1 ive cryogenic 
tankage. Tliis should be done to provide a basis fur comparison 
with similar tests which have already been conduetod by the I'A/V 
nsing surplus P 2 alrcvuit with eouventional fuels. 

5. Demonstrate the e£feeLivenu.ss of purge systems to eliminate 1 1 re 

hazard in even (; of onboard fuel leaks, especially witli the cryogenic 
tuuls. Also, demonstrate the ef fectivenesu of candidate systems to 
detect the presence of fuel leaks. 
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AlTKNPiX A 


hmirncM, uiitiUj.TH oi'' xiU'; aihcuai-x uAPiiOUti 

l)itJPKU!iJ.ON MOIIKJ. 

uucli t)l' tlu> licuiiuriotJ iXuLcul jLu tabXu 2*), ihuj'e ia a ph>i. oi 
uiLitudu (Z) vurauu diauaucu (X) »1 thu U):aji:ctoi.’y ul tlic ,larj.'‘-;«t- pall , Tlio 
timu ill uuuuiidu uu which ihc ^aacoua i'ucl cuiicuiitraLicm drupa below iui JowcoL 
liluimiiublc liiiiil. ia plotted at the end oi the trajectory. Each plot aliowa the 
iucl, ita lowuat i'luiiuiiuble mole iractioa , circular apill radluti or axial 
upill width ill meturut wind apeed at 10 metura height la metcra per Hecuiidt 
auri’acc temperature iu kclviu, and relative humidity (dimeiiaioiileHw) . 

Each plot ia tollowed by a tubular liatiiig of time, dc.wawlud diataiice 
of pul'i center I pul'l' temperature, atmoapheric temperature at pui'i height miiiub 
pui't' temperature, buoyancy acceleration, pul’l: vertical apeod, pull: lieight above 
ground, overall angle of aacctit of puff from aturting point, puff volume, 

3 uutrainment speed terms, gaseous fuel mole fraction in puff, incremental 
entrainment volumes through puff sides, tup and bottom, radius or circular 
puff, half thickness of puff, and the length and width of an axial puff. 

The listing also contains the input conditions: the circular spill 

radius or the axial spill width in meters, the fuel, the surface elevation 
in meters above sea level, the evaporation rate in meters per second, the 
surface temperature in kelvin, the relative humidity, the atmospheric 
stability class number, the wind speed at 10 meters above the ground in meters 
per second, and the gaseous fuel lower flaiiuiiablu coueentratlor expressed as 
u mule fraction. 

All scenarios in this report were run at zero surface elevation (sea 
levels and surface temperature of 294®K. Also, all scenarios except for 
extra runs of Scenario 4A with.Lll2 were run with the most stable atmosphere 
(class b) and 50% relative humidity ( «• O.S). 
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TAUI.K .'.U. ■ tiAtiKUUli l)J.‘il'i;U;iiUN JiCKNAHitW AMAIA7.1';)) 


tiiiil 

Maxliimm 

Kadliiii 

(m) 

Maxliiimn 
AkIuI Width 
(nil 

hvaimratidii 

Hate 

liii/i) 

III;, 

u;ii^ 

1 

1,3 


2.24 K 10^ 
U.Ofi X to'’ 

Jl''4 

iti.b 

•- 

2.00 X 1 U'*® 

Jet A 

utt 


7.00 X 10 “ 

Ul;, 
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2 .b 0 x 10 '^ 

LCII 4 

3,4 


4.00 X 1 o'* 

JIM 

40 

~ 

2.00 X 10 '** 

Jut A 

212 
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JP4 
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Jut A 
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LH 2 
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Jut A 
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LII 2 
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LGH 4 
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JP-4 
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Jet A 
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LH 2 
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LCII 4 
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Jet A 
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LMj 
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ICH 4 
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1.62 X lO ** 

JP-4 

74 
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2 . 00 x 10 ® 

Jet A 

171 


7.00 X 10® 


(Joiiimiiiitii 


Nu fliimmoblu cuiicttiitmtlun 


Nu (lamiMBljlu cu'icuntratiun 


Nu flemmalilu cuiicuiitratiun 


Nu flammablo curiceiit'ation 


No tlemmablu uuticuiitratiutt 


Nu flaminable vuDcaiitratiun 


Nu flaiMinablu cuncautratiUM 


Nu flammable cunuaiitratiuii 



Ill .ulUitJ.im to Uu'. alinviii l-hu l.iiat. twu aurimi oJ' aompui.o: rt*,fji 4 ii.a 
111 ApiHMidlx A arn loi: i.lu> rn,l l.awliiy iiaiieii liivu,l.vJ.iiy ficuiiar:J.ii -'lA w.Uh Ul, 
i. iiu i : 

• h.lx i-UKK.iii , all his 'll t. J 4^,1 1 , tixtiupi; foi* ViiKiaUiiimi ol atiiiofipluirlc; 
fUiibi l .l.ty . Till, riiiiri ui-v .lluced iu ordur of cXahh mi tolliiwas 0, !i , 

4, ;i, Z, mid .1., 

• l'’Xv« ca.'imi, a], I liluiitlua.l., exucipt foi' variationa oJ. relative huiiiidity, 
Tliu niiia apprar lu order; 0. lio, 75 and 100 pereent. 
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